arXiv:1508.04869v2 [quant-ph] 24 Mar 2016 


Ultrafast Optimal Sideband Cooling under Non-Markovian Evolution 


Johan F. Triana, Andres F. Estrada, and Leonardo A. Pachon 
Grupo de Fisica Atomica y Molecular, Instituto de Fisica, Facultad de Ciencias Exactas y Naturales, 
Universidad de Antioquia UdeA; Calle 70 No. 52-21, Medellin, Colombia. 

(Dated: March 25, 2016) 

A sideband cooling strategy that incorporates (i) the dynamics induced by structured (non- 
Markovian) environments in the target and auxiliary systems and (ii) the optimally-time-modulated 
interaction between them is developed. For the context of cavity optomechanics, when non-Markovian 
dynamics are considered in the target system, ground state cooling is reached at much faster rates and 
at much lower phonon occupation number than previously reported. In contrast to similar current 
strategies, ground state cooling is reached here for coupling-strength rates that are experimentally 
accesible for the state-of-the-art implementations. After the ultrafast optimal-ground-state-cooling 
protocol is accomplished, an additional optimal control strategy is considered to maintain the 
phonon number as closer as possible to the one obtained in the cooling procedure. Contrary to the 
conventional expectation, when non-Markovian dynamics are considered in the auxiliary system, the 
efficiency of the cooling protocol is undermined. 
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Introduction .—Fabricating and controlling micro- and 
nano-devices is of primary importance to develop quantum 
technologies. In this effort, optomechanical devices play 
a central role because they may become a building 
block for hybrid systems [1-4]. To benefit from 
quantum effects, potential quantum technologies that 
utilize optomechanical systems need to prepare the 
mechanical component in the ground state. Thus, 
searching alternative schemes and refining present 
strategies to improve cooling protocols are of great interest. 
However, the highly non-trivial interplay between time- 
dependent external fields and structured environments 
poses constraints on the coherent control of these systems 
[5] that have been barely explored. 

Despite the copious recent literature on optomechanical 
cooling [1-3, 6-12], a key point of the physics at 
low temperature remains unexplored, namely, low 
temperature induces non-Markovian dynamics [13-16]. 
This very fact and the recent experimental evidence that 
the dynamics of micro-resonators are non-Markovian 
[17], certainly, suggest that the present understating 
of cooling processes is incomplete. Moreover, because 
these approaches work on very short time-scales, of 
the order of the period of the resonator, it is expected 
that non-Markovian effects, either from low-temperature 
fluctuations or structured environments, dominate the 
energy and entropy transfer. 

The exploration of non-Markovian dynamics has 
already leaded to the enhancement of, e.g., quantum 
speed limits [18], survival of entanglement [16, 19, 20], 
improvements in quantum metrology [21] and corrections 
to thermal equilibrium states [15, 22]. In the context 
of cooling, it was even shown that only under non- 
Markovian dynamics the entropy of a parametrically 
driven resonator can decrease [5]. The non-Markovian 
aspects of the dynamics discussed above are introduced 
here to one of the most promising techniques to reach the 


minimum phonon number, namely, finding an optimal 
coupling function via optimal control theory [5, 8]. Based 
on analytic exact results derived in the context of the 
Feynman-Vernon influence theory, it is found that when 
the non-Markovian character is considered in the target 
resonator, the phonon number is lower than the predicted 
by Markovian processes; however, in contrast to the 
conventional expectation, non-Markovian dynamics in 
the auxiliary system deteriorates the cooling protocol. 

Model and Theory .—For mechanical systems, sideband 
cooling consists in coupling a resonator to a microwave 
or optical resonator whose frequency is sufficiently high 
that it effectively sits in its ground state at ambient 
temperature. The coupling between the mechanical 
and electromagnetic modes is mediated by radiation 
pressure and is nonlinear. However, under realistic 
experimental conditions [9, 23-25], the coupling can be 
assumed to be linear [2, 4]. This common model with 
linear coupling is considered below [? ]. Specifically, 
in the linear approximation, the system consists in two 
harmonic modes that describe the mechanical mode and 
the optical mode with the annhiqulation operators d and 
b, respectively. The Hamiltonian reads 

H = HuJiTirnO^d + Hujonib^b + flg{t){d^ -j- a)(6’^ -b b) (1) 

with Wmm and Wom, the frequency of the mechanical 
and electromagnetic modes, respectively, and g{t) is 
the arbitrary time-dependent optomechanical coupling 
function to be found via optimal control theory [26]. For 
later convenience, the period of the mechanical mode is 
labeled by rmm = 27r/a;mm- To describe the interaction 
with the environment of the resonator and the losses in the 
cavity, each mode is coupled to an independent thermal 
bath with arbitrary spectrum and described in the context 
of the Ullersma-Caldeira-Leggett model [27, 28] [see the 
Supplementary Material for details]. 
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The particular functional form of the Hamiltonian (1), 
the environment model and the thermal initial states 
allow for the complete description of the dynamics in 
terms of the variances of the modes’ coordinates (see, 
e.g., Ref. [16]). This was utilized to combine sideband 
cooling and optimal control theory [8]. Specifically, the 
dynamics were formulated in terms of an adjoint master 
equation for the quadratures of the modes derived from 
the Markovian version of the Brownian-motion master 
equation [29, 30]. However, this master equation may 
violate the density-operator positivity (see Sec. 3.6 in 
Ref. [31]). Moreover, for time dependent Hamiltonians, 
an adjoint Lindblad master equation can be derived only 
if the Liouvillian of the non-unitary dynamics commutes 
with its associated-time-ordered propagator (see Sec. 3.2 
in Ref. [31]), this is not the case for the Hamiltonian (1). 

To circumvent the issues raised above, the Feynman- 
Vernon’s influence functional theory [32] is employed here. 
It allows for performing a description of the dynamics 
without any approximation [see Supplementary Material], 
and it further obeys the positivity of the density operator 
and straightforwardly incorporates driving fields. By 
combining the exact analytic results provided by the 
influence functional approach and an efficient numeral 
algorithm to solve integro-differential equations, the 
dynamics are solved here for arbitrary bath spectrum 
and driving field. 

Non-Markovian Optimal Sideband Cooling .—The po¬ 
tential benefits from cooling optomechanical systems have 
attracted a great deal of attention [1, 2, 7, 10, 23, 33, 34] 
and recently, it has been powered by optimal control the¬ 
ory [3, 8]. However, a key dynamical feature has been left 
out of the discussion, namely, the non-Markovian dynam¬ 
ics induced by low temperature and structured environ¬ 
ments [13-16]. Therefore, the goal here is to calculate and 
investigate the dependence of the minimum phonon num¬ 
ber in the mechanical resonator (h)(tcooi), at the shortest 
possible time tcoop on the non-Markovian character of 
the dynamics. In terms of the second moments of the 
position [q^) and momentum (p^), the phonon number at 
time t is given by 




(P")(i) 


m 


rnumm {q )it) 



The optimization procedure under non-Markovian dynam¬ 
ics, described in the Supplementary Material, leads to 
the minimum number phonon occupation at time tcooi- 
The cooling time tcooi is chosen as short as possible by 
hand. The spectrum of the mechanical-resonator’s envi¬ 
ronment is described by the spectral density Jmm(w) = 
7 Wq/ (wq + whereas the cavity environment is cho¬ 
sen, for convenience, as Jom(w) = kIIq/ (Hq As 

in Ref. [8], it is assumed that the frequency conversion 
is exact and set Wmm = Wom- The corrections to this ap¬ 
proximation are of the order of (wmm/wom)^ This spectral 



FIG. 1. (a) Time evolution of the phonon number n{t) during 
the optimal modulation of the interaction with ric = 0 , 7 = 
10“®cum and ft = 2.15 x lO^^cum. In particular, for nx = 10'^ 
(nx = 10 ®) and final time tcooi = O.SSxm (Gooi = 1.25Tm), the 
dashed red (blue) curve depicts the dynamics of the phonon 
number under Markovian dynamics (cud = Hd = lOOtUm) while 
the continues black (magenta) curve depicts the dynamics 
under non-Markovian dynamics (tUD = Hd = tUm). The 
continues dark green (cyan) curve depicts the non-Markovian 
phonon-occupation-number dynamics for nx = 10 ® (nx = 10 ®) 
and tcooi = 1.25rni (tcooi = 0.55rni) whereas the dashed olive 
(light green) curve stands for the Markovian dynamics. The 
y-axis is in logarithmic scale, (b) Optimal coupling function 
with tcooi = 0.55rni under non-Markovian (continuous lines) 
and Markovian dynamics (dashed line), (c) Optimal coupling 
function with Gooi = 1.25Tm under non-Markovian (continuous 
lines) and Markovian dynamics (dashed line). 


density relates to temperature fluctuations [? ]. 

For typical values of optomechanical setups, Fig. 1(a) 
shows the optimal dynamics of the phonon number in 
the Markovian (dashed lines) and in non-Markovian 
(continuous lines) regime. Figure 1 shows two alternatives 
to decrease the minimum phonon number, namely, (i) 
extending the cooling time or (ii) considering non- 
Markovian effects in the dynamics. By lengthening the 
cooling process from tcooi = 0.55Tm (black and cyan 
lines) to tcooi = l-25Tm (green and magenta lines) and 
for the initials phonon number in the mechanical mode 
riT = 10® and nr = 10®, it is possible to decrease the 
minimum phonon number by an order of magnitude, 
from ca. (n(tcooi)) = 10“® to (n(tcooi)) = 10“®. The 
Markovian and non-Markovian character of the dynamics 
does not change considerably the functional form of the 
optomechancal coupling; however, the small difference 
is enough to decrease the minimum phonon-number by 
around one order of magnitude (see also Table I). 

The coupling strength for the ultrafast cooling 















3 



(7z(tcool)) 

7/wm 

opt. omM+ mmM 

om„M+ nim„M 

opt. omnM+ mmnM 

omM+ mmnM 

opt. omM+ mmnM 

lO-*^ 

9.03 X 10-^ 

8.86x10-3 

3.43x10-3 

1.91x10-3 

9.99x10-3 

10-'^ 

1.04 X 10"^ 

1.02x10"^ 

4.79x10-3 

2.60x10-3 

7.91x10-4 

10”4 

3.28 X 10"^ 

2.40x10"^ 

1.99x10-^ 

1.63x10-^ 

1.52x10-^ 

10-^ 

2.61 X 10-1 

1.61x10-3 

1.53x10-3 

1.53x10-3 

1.51x10-3 

1 

0 

1—1 

2.45 

1.52 

1.50 

1.50 

1.50 

10-^ 

21.12 

14.21 

13.52 

13.52 

13.52 


TABLE I. Minimum phonon number in the mechanical resonator at a time tcooi = 0.55rni for different values of the dissipation 
rate 7 /aJm and for different cooling scenarios (see text). Cavity dissipation rate k = l0~*LOm- 


scenario in Fig. 1(b) is beyond the present experimental 
capabilities; however, results for the fast cooling case in 
Fig. 1(c) are encouraging [23]. Moreover, in the regime of 
Fig. 1(c) cooling from room temperature is possible [12]. 

The non-trivial interplay between non-Markovian 
dynamics and driving fields is explored in Table I, it 
displays the predicted phonon number for a variety of 
scenarios, namely, (i: opt. omM+ mmM) the minimum 
phonon number obtained from optimization under 
Markovian dynamics, (ii; om„iy[+ mninM) the phonon 
number obtained under non-Markovian dynamics with 
the optimal coupling found under Markovian dynamics, 
(iii; opt. omnM+ nimnM) the phonon number obtained 
from optimization under non-Markovian dynamics, 
(iv: omM+ mm^iu) the phonon number obtained under 
Markovian dynamics in the optical mode and non- 
Markovian dynamics in the mechanical mode with the 
optimal coupling found under Markovian dynamics and 
(v: opt. omM+ nim„M) the phonon number obtained from 
optimization under Markovian dynamics in the optical 
mode and non-Markovian dynamics in the mechanical 
mode. 

Comparison of scenarios (i) and (ii) in Table I shows 
that, for the same coupling function, the presence of 
non-Markovian dynamics reduces the minimum phonon 
number and that this decrease is more noticeable when the 
dissipation rate of the mechanical resonator increases. The 
third scenario depicts the influence of the optimization 
process. Even though there is a reduction in the phonon 
number respect to the second scenario, in absolute terms, 
the decrease is tiny. Besides, if the optical-mode dynamics 
are considered as Markovian and the dynamics of the 
mechanical mode as non-Markovian, fourth scenario, there 
is a decrease in the number of phonons for low values of the 
decay factor y/wm- Surprisingly, only when Markovian, 
non-Markovian dynamics and optimal control theory are 
combined, fifth scenario, a substantial reduction in the 
phonon number is reached for low mechanical decay rates. 
For these low-decay-factor cases, which are relevant under 
experimental conditions, the strategy in the fifth scenario 
is able to bring the number of phonons more than one 
order of magnitud below the number of phonons reached 
under the conditions of the first scenario. 


This non-trivial and unexpected interplay between 
Markovian, non-Markovian dynamics and optimal control 
theory is also observed for a variety of values of the cavity 
decay rate in Table II. The advantage of a Markovian- 
dynamics cavity is more transparent for large values of 
the cavity losses rate k, e.g., for a non-Markovian cavity 
(Hd = Wc) with K = 2.15 X lO^^Wc, there is no ground- 
state cooling [(unM(tcooi) = 2.34]. 


^cool / 

k/Wc 

('^M (^cool)) 

(^nM(^cool)) 

('n-OIIlM (^cool)) 

0.55 

1x10-3 

0.016 

0.015 

0.014 

0.6 

1.5x10-2 

0.025 

0.019 

0.019 

0.8 

2.5x10-2 

0.029 

0.030 

0.021 

0.8 

4.5x10-2 

0.035 

0.060 

0.026 

0.8 

5.5x10-2 

0.037 

0.086 

0.032 

1.0 

1.25x10-3 

0.048 

0.356 

0.040 

1.6 

2.15x10-3 

0.056 

2.34 

0.044 


TABLE II. Minimum phonon number for scenarios 
(i: opt. omM-f mmM) (hM(tcooi), (hi: opt. omM+ mmM) 
(hnM(tcooi) and (v: opt. omM+ mmM) (uomM(tcooi) for dif¬ 
ferent cooling times and for different parameters of cavity 
dissipation. The initial parameters are nx = 100, 7 = 


Because the mechanical-mode initial state is highly 
thermally-populated, the non-Markovian character of 
the dynamics originates mainly from the structure of 
the thermal bath and not from quantum fluctuations at 
low temperature. This explains why, for non-structured 
environments, Markovian descriptions of the dynamics 
may had provided sensible results. 

Maintaining the minimum phonon number over 
time .—Reaching a very low phonon number in a short 
period of time is a desirable goal. However, because the 
resonator is continuously coupled to its environment, keep¬ 
ing that phonon number is a must. In doing so, introduce 
a second optimal control scenario where the objetive is 
to maintain, over a long period of time, the minimum 
number of phonons in the resonator reached in the ultra¬ 
fast cooling in Fig. 1. Define then a composite optimal 
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FIG. 2. (a) Phonon number as a function of time during the 
optimal control protocol aimed at maintaing the minimum 
phonon number for different values of the dissipation rate. The 
initial parameters are nx = (u(tcooi)), ric = 0, k = 2.15 x 10“^. 
(b) The optimal optomechanical coupling function c{t) for 
7 = (c) The optimal optomechanical coupling 

function gni{t)- 


optomechanical coupling function as 


f 5c(i) 0 < t < tcooi, 

\ ^ ^ ^cool; 


( 3 ) 


which encompasses the entire cooling process, i.e., 
reaching a minimum phonon number by applying gdt) and 
maintaining this number once it is obtained with 
Because performing the second optimization process 
under non-Markovian dynamics over long times requires 
an enormous amount of computing time, the second 
optimization process is performed under the assumption 
that for long times, the Markovian approximation holds 
so that the optimal coupling function is calculated for 
Markovian dynamics. 

Figure 2(a) depicts the time evolution of the phonon 
number under the action of the optimal coupling function 
c{t) for a variety of experimentally-relevant values of the 
decay rate 7 . To obtain an experimentally-accessible- 
coupling-strength scenario, |c(t)| < 10 “^, the cooling 
time is set to tcooi = 6 rm and the initial phonon number 
to riT = 10^ [see Fig. 2(b)]. For the second control phase, 
the initial number of phonons is set to be the minimum 
phonon number reached at tcooi- To maintaining the 
minimum phonon number requires moving from the strong 
coupling regime to a weak coupling regime between the 
two modes, see Figs. 2(b) and 2(c). Likewise, with the 
optimal optomechanical coupling function c{t) the phonon 
number is approximately maintained for 50 periods of the 
mechanical resonator. 

Discussion. —Analytical solutions and an efficient op¬ 
timal control protocol showed that the presence of non- 
Markovian dynamics allows for lower phonon numbers 
than the predicted by Markovian-and-KWA-based pre¬ 


vious works. Surprisingly, significant enhancements are 
found as the interplay between non-Markovian dynamics 
in the mechanical mode, Markovian dynamics in opti¬ 
cal mode and optimally-designed coupling functions. To 
understand this effect, note that for well-behaved J(w), 
non-Markovian dynamics define an effective coupling to 
the thermal bath [15, 16] that, in general, is weaker than 
the Markovian one. Thus, when non-Markovian dynamics 
are considered in the cavity dynamics, the rate at which 
the electromagnetic mode releases the entropy, that ab- 
sorbes from the mechanical mode, into its environment 
diminishes and therefore, the number of phonons in the 
mechanical mode does not largely decrease compared to 
the case of Markovian-dynamics cavity. Furthermore, be¬ 
cause the mechanical mode is constantly coupled to its 
environment, the optimal cooling protocol was extended 
[see Eq. (3)] to preserve the phonon number as low as 
possible after cooling is reached. For this second protocol, 
the coupling amplitude is between reach of present tech¬ 
nology and can be combined with already-experimentally- 
implemented cooling protocols [9, 23-25]. 

The optimization protocol developed here relies on the 
optimization over the Green functions of the trajectories 
that minimize the influence functional (see Supplementary 
Material) and can be readily implemented in semiclassical 
formulations of quantum mechanics in phase space [35, 36]. 
Thus, non-linear systems can be addressed and the 
influence of quantum fluctuations in the design of optimal 
pulses of coupling functions can be analyzed. This enables 
the present proposal, e.g., in the context of optimal control 
theory of molecular processes [37]. 
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INFLUENCE FUNCTIONAL 


The state of the mechanical mode and the optical mode q2+, Qi-, (12,^) can be 

determined from 


p{q”+, q2+, q”-, q2, t)= dq[^dq2+dq[_dq2_J{q”+, ^ 2 +, q”-, q'L^t'^ 0 ) 


X ps{q'i+,q'2+^q'i-,q'2-^d) 


( 1 ) 


where ps('?'i+, 92 +^ ?2-> 0 ) accounts for the initial state. J (g"+, ^2+, q'L, q'L^t'^ ?'i+> ^2+) ?2-> 0 ) 

is the propagating function to be obtained from the Feynman-Vernon influence functional ap¬ 
proach. After assuming that each mode is coupled to its own Ullersma-Caldeira-Leggett-type 
bath at temperature Ti^2, the propagating function reads 


7 / //////// 4 . / / / / r»\ 

J (gi+, q2+, qi-: 92-, ?2+, Ql-: ^2-, 0) = 


1 r f'p 

—— exp < - ds 

N{t) [hJo 


X exp < — - / ds 

^ Jo 


0=1 

0 = 1 


- c(s)gi+(s)g 2 +(s^ 
- c(s)gi_(s)g 2 -(s) 


X J^[qi+,q2+,qi-,q2- 


( 2 ) 


with 


2 

J^[qi+, q2+, qi-,q2-] = n ®xp 

0=1 


im, 


2 ^ I (?«+ + ds7„(s)[g„+(s) - go-(s)] 


+ ds du'yds - u)[qa+{u) + qa-{u)][qa+{s) - qa-{u)] 

Jo Jo 

xexpj-i/ ds [ duKa{u - s)[qa+{s) - qa-{s)][qa+{u) - qa-{u)]} , 


h 


'0 Jo 


( 3 ) 


2 



where 


7a (s) 

Ko.{s) 


dUoi Jai^^a) 


m 


a Jo 

“du;. 


vr Uln 


cos (UaS), 


TT 


Ja(a;a)coth ( ] cos(a;Qs), 


( 4 ) 

( 5 ) 


being Ja{oJa) the spectral density. In the main text, Ja{oJa) was chosen such that is accounts 
for Ohmic dissipation with a cutoff frequency uJod- The functions qa±{s) obey the following 
relations (1 for mechanical mode and 2 for optical mode) 


/ N 2 / N c(s) , , Id 

5 'l+(' 5 ) +'^lQ'l+('S) H- (l2+{s) + - — 

mi 2 ds 


f't pt 

dM7i(s - M)(gi++ gi_) - / dM7i(M - s)(gi+- gi_) 


c(s) 


Qi-is) +ujiqi-{s) H- q 2 -is) + 77- 

mi 2 ds 


L-'o 

1 d r /-i 


I duji(s-u)(qi+ + qi_)+ duji(u-s)(qi+-qi_) 

0 Js 

q2+(s) + ufq2+(s) + ^gi+(s) + ^ dM72(s - u)(q2+ + q2-) - du'y2{u - s)(g2+ - ^2-) 

i dM 72 (s - M)(g2+ + ^2-) + / dK 72 (M - s)(g2+- g2-) 

0 Js 

( 6 ) 


m 2 


\ 2 / \ c(s) , , Id 

Q2-is) + uj;^^q2-{s) H- qi-is) + - — 

m2 2 ds 


= 0 , 
= 0 , 
= 0 , 
= 0 . 


To simplify the above expressions and subsequent calculations, introduce the half sum 
and difference coordinates as 


Qa 2 ^^®+ 4“ Qa—')i Qa qa+ qa— 


( 7 ) 


where the Jacobian of the coordinates transformation is equal to one. 

The interest is to hnd the minimum phonon number in the mechanical mode, which in 
terms of the second moments is given by 


{n){t) = 


2fkv 


{p^m 


m 


+ mu {q ){t) 


1 

2 ' 


( 8 ) 


Following the technical details in Ref. [ 1 ] , the second moments read 


(Qi) (t) = 

(pI) (t) = I dQi Qu t) 


dg;'g;'pi(g;',g'; = o,t), 

roo ^2 


dq'l 


g "=0 


( 9 ) 

( 10 ) 
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with 


Pi(gi+,gi_,t) = {qi+\tT 2 p{t)\q^_) = / dQ2p{qi+,Q2,qi-,q2 = 0,t). 


( 11 ) 


EQUATIONS OF MOTION AND ITS SOLUTION 


The differential equations of motion for the coordinates Qi and qi take the form 

Qi,2{s) + ul2Ql,2{s) ^ -1Q2,i(s) + — / dw 7i,2(s - m)Qi,2(m) = 0, 

mi,2 ds Jo 

2 c(s) d /”* 

gi,2(s)+a;i2gi,2(s) H-g2,i(s) - — / dn7i,2(M - s)gi,2(M) = 0 , 


mi,2 


and have the following boundary conditions 


ds 


. . Q', s = 0 I g', s = 0 

Qi[s) = { , qjs) = \ 

Qi, s = t q'l, s = t 


(13) 


The solution to Eqs. (12) can be expresed as 


Qi(Jy 5) — Ui(t, s)Qi + U2(t, s)Q” + U^ij, s)Q 2 + U^ij, s)Q 2 
Q2{t, s) = Vi{t, s)Q2 + V2{t, s)Q2 + V3(t, s)Qi + V 4 (t, s)Q'l, 
qi{t, s) = Ui{t, s)q[ + U2{t, s)q'l + Usit, s)q'2 + Ui{t, 5)^2, 
q2{t,s) = Vi(t,s)q'2 + V2{t,s)q2 + Vo{t,s)(J^ +v^it, s)q'{, 


(14) 


where the boundary conditions that satisfied the functions Ui, Vi, Ui and Vi in Eqs. (14) are 
given by 

Ui(t,0) = l, Ui(t,t) = 0, U2it,0) = 0, U2(t,t) = l, Ui{t,0) = Ui{t,t) = 0, (15) 

Ui(t,0) = l, Vi{t,t) = 0, U2(t,0) = 0, V2{t,t) = l, Vi{t,0) = Vi{t,t) = 0, (16) 

Ml(t,0) = l, Ui{t,t)=0, U2(t,0) = 0, U2it,t) = l, Ui(t,0) = Ui(t,t) = 0, (17) 

ni(t,0) = l, vi{t,t)=0, V2{t,0) = 0, V 2 {t,t) = l, Vi{t,0) = Vi{t,t) = 0, (18) 


for i = 3,4. 
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One form to solve the problem is solving in terms of the initial conditions 

Q*(0) = g*(0) = 1, 4(0) = g,(0) = 1, g,(0) = g,(0) = 0, and 4(0) = g,(0) = 0, 

where Ui(t,s), Vi{t,s), Ui(t,s), Vi(t,s) in terms of the solutions (pi{s), (f>i{s), r'j(<s) and ^i{s) 
(solution functions of the problem initial condition) are given by 


Ui{t,s) = ipi{s) - 


+ 


<^l(f)V? 2 (s) 02 (t) 


9^2 (t) 03(f) 9 ^ 4 ( 5 ) 


9 ^ 2 (f) 02 (f) - 9^4(f)04(f) V^ 2 (f) 02 (f) “ V^4(f)04(f) 

V?2(g)03(f)<^4(f) V?i(t)9?4(-S) 04(f) 

V^ 2 (f) 02 (f) - V^4(f)04(f) V^ 2 (f) 02 (f) “ 9?4(f)04(f) ’ 


(19) 


U2{t,s) = 


<^ 2 (s) 02 (f) 


<^4(-S)04(f) 


V^ 2 (f) 02 (f) - V^4(f)04(f) V^ 2 (f) 02 (f) “ V54(f)04(f) ’ 


( 20 ) 


U3{t,s) = 993 ( 5 ) 


<^2(s)02(f)<^3(f) 


01 (f) 9^2 (f) <^ 4 ( 5 ) 


+ 


9^2(f)02(f) - <^4(f)04(f) (P2{t)h{t) - 9^4(f)04(f) 

0l(f)V?2(-S)9^4(f) V?3 (f) 9^4 (g) 04(f) 

V^2(f)02(f) - V94(f)04(f) V^2(f)02(f) “ 994(f)04(f) ’ 


( 21 ) 


Ui{t,s) = 




<^2(s)04(f) 


9^2(f)02(f) - V^4(f)04(f) V92(f)02(f) “ V^4(f)04(f) ’ 


( 22 ) 


l/l(f,s) = 01 ( 5 ) - 


+ 


0 l(f)<^ 2 (f) 02 (s) 


02(f)<^3(f)04(-S) 


<^2(f)02(f) - 994(f)04(f) V^2(f)02(f) “ V94(f)04(f) 

0l(f)994(f)04(g) 02(g) 973 (f) 04(f) 

V72(f)02(f) - 974(f)04(f) 972(f)02(f) “ 974(f)04(f) ’ 


(23) 


V^ 2 (f,g) = 


<d 2 (f) 02 (g) 


974(f)04(g) 


972(f)02(f) - 974(f)04(f) 972(f)02(f) “ 974(f)04(f) ’ 


(24) 
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V3{t,s) = 03(s) - 


+ 




<^lW02(t)04(-S) 


<^2(t)02(t) - <^4(t)04(t) ip2it)(t)2{t) - 994(t)04(t) 

03 (t)v?4(t) 04(g) 9?l(t)02(g)04(t) 

V^ 2 (t) 02 (t) - 924(t)04(t) V^ 2 (t) 02 (t) - (p ’ 


( 25 ) 


I4(t,g) = 


V? 2 (t) 02 (g) 


(pA{t)(t)i{s) 


<P2{t)4>2{t) - <P4{t)4>4{t) P>2{t)4>2{t) - (pA{t)(j)4{t) ’ 


(26) 


Ui{t,s) = Ui{s) 


Ul{t)u2{s)^2{t) 


l^2{t)^3{t)Ms) 


^^2(g)^3(t)z24(t) Z/i(t)z/4(g)^4(t) 

jy2{t)'&2{t) - iy4{t)'&4{t) Jy2{t)'&2{t) - iy4{t)'&4{t) ’ 


(27) 


U2(t,s) = 


U2{s)^2{t) 


iy4{s)^4{t) 


T^2it)'&2it) - V4it)'&4it) U2{t)'&2it) - V4it)'&4it) ’ 


(28) 


U3{t,s) = Z/ 3 (s) - 


+ 


V2{s)^2{t)v3{t) 


^l{t)v2{t)v4{s) 


T^2it)^2it) - P4it)^iit) T^2it)^2it) - P4it)^4it) 
^l{t)v2{s)v4{t) Z/3(t)z/4(g)^4(t) 

l^2{t)^2{t) - V4{t)^4{t) iy2{t)'&2{t) - iy4{t)l!^4{t) ’ 


(29) 


U4(t,s) = 


l^ 2 {t) 124 ( 3 ) 


I22{s)^4{t) 


^2{t)'&2{t) - 124{t)'&4{t) ^2{t)'&2{t) “ 124{t)'&4{t) ’ 


(30) 


Vi{t,s) = ^i{s) 


^l{t)l22{t){}2{s) 


^2(t)l23(t){}4{s) 


T22{t)d2{t) - l'4{t)d4{t) I22{t)d2{t) - l'4{t)d4{t) 

^l(t)z/4(t)^4(g) ^2(g)^^3(^)^4(t) 

T22{t)d2{t) - l'4{t)d4{t) l'2{t)d2{t) - l'4{t)d4{t) ’ 


(31) 


V2{t,s) = 


V2{t)^2{s) 


V4{t)^A{s) 


T^2it)'&2it) - U4{t)'&4{t) U2{t)'&2it) - V4it)'&4it) ’ 


(32) 
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(f ^ ,Q _ t/i(t)^2(t)^4(g) 

^ ^ Mt)^2{t) - iy4{t)^dt) Mt)^2{t) - iyA{t)i^4{t) 

l^2it)'&2{t) - I24{t)l^4{t) l^2it)l^2it) - I24{t)l^4{t) ’ 


V4it,s) 


I22{t)-d2{s) _ U4it)^4{s) 

V2{t)d2{t) - V4{t)d4{t) h'2{t)d2{t) - V4{t)d4{t) ' 


( 34 ) 


Therefore, in terms of the solutions of the initial condition problem with ipi, (f)i, Ui and 'dt, 
the equations of motion are transformed as 


cftsl d. 

<fi,3{s) + ulipi^3{s) H-^(/>3,i(s) + q- / du7i(s - u)ipi^3{u) = 0, 

mi ds Jo 

/ \ 2 / N c(s) , / . d , , , , , 

^P2,4{S) + UJ^^P2,4{S) H-04,2('S) + U” / dM7i(s - U)^p2,4{u) = 0, 

mi ds 


^2,4 


'0 
d 


(35) 


3 ( 5 ) + 02201, 3 ( 5 ) H-^‘^3,l(s) + 7 - / dM72(s - m)0i,3(m) = 0, 

m 2 ds Jo 

cfsl ci 

,4(s) + 02202, 4 ( 5 ) H- -^4,2{s) + 7 - / dM72(s - m)02,4(m) = 0, 

m 2 ds Jo 


hl,3(s) +02 iZ/i,3(s 

i^2,4(s) + t(2^Z/2,4(s 
i9i,3(s) + 022^9i,3(s 
'^ 2 , 4 ( 5 ) + 022'd2,4(’S 


c(g) 

mi 

c(5) 

mi 

c(g) 

m 2 

c(g) 

m 2 


■^ 3 , 1 ( 5 ) 
^4,2(s) 
^3,1 (s) 
^ 4 , 2 ( 5 ) 


d 

ds 

d 

ds 

ds 

d 

ds 


dti 7i (m 


du 7i (m 


dM72(M 


dM72(M 


s)i/i, 3 (m) = 0 , 

s)u2,4{u) = 0 , 

s)'d2,4(M) = 0 . 


(36) 


In this work, the dissipation kernels in Eqs. (35-36) for the mechanical and electromagnetic 
mode are given by 

71 ( 5 ) = 7l'^lDe■'^l°^ 72 (s) = 7202206 "^^™*. (37) 

The equations of motion to solve and optimize are obtained after replacing these dissipation 
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kernels in Eqs. (35-36) and are given by 


01,3(5) -|-02^(^i,3(s) -|- 

S 0 

CO 

-U 

- 01,2(5) — 0, 

01,3(5) -1- 02201,3(5) -f 

C(5) . . 

9^3,1 5 d 

m2 

- 03,4(5) = 0, 


ji(s) -t- a;iDJi(s) - 7ia;iDV7i(s) = 0 , (' 33 ) 

'^ 2 ( 5 ) + WiD ^ 2 ( 5 ) — 71^10(^3 (s) = 0, 

^ 3 ( 5 ) 0220 ^ 3 ( 5 ) - 72i:^2D0i(s) = 0 , 

^ 4 ( 5 ) + (^2YiJi{s) - 72022003(5) = 0, 


02 , 4 ( 5 ) -\- 0j‘l^2,i{s) -\- 


C(5) 

mi 


5 , 6 ( 5 ) — 0, 


2 c(s) • 

02 , 4 ( 5 ) -|-02202,4(5) H-(^ 4 , 2 ( 5 ) -|- Jlfiis) = 0, 

m 2 

j^^s) + 0210 ^ 5 ( 5 ) - 71^1D</?2(5) = 0, 
06(5) -+• 0210 ^ 6 ( 5 ) - 71 ^ 109 ^ 4 ( 5 ) = 0, 
^ 7 ( 5 ) + 022D</7(5) — 72022002 ( 5 ) = 0, 
08(5) 0220 ^ 8 ( 5 ) - 72022004(5) = 0, 


(39) 


c( 5 ) 

^^ 1 , 3 ( 5 ) -f 02 ^ 1^1 3 ( 5 ) H-'*^ 3 , 1 ( 5 ) —^’ 1 , 2 ( 5 ) = 0 , 

mi 

'*^1,3(5) -f 022191,3(5) -|- ~~^ 3 ,i( 5 ) —73,4(5) = 0 , 

ji{s) - Uiojiis) + 7 i 02 idi ^ i ( 5 ) = 0 , ^ 40 ) 

72(5) - 021072(5) -f 7 i'^io 1^3(5) = 0 , 

73(5) - 022073(5) + 72022 Oi 9 i( 5 ) = 0 , 

74(5) - 022074(5) + 720220193(5) = 0, 
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(41) 


with 


c( 

i>2,4(s) +ojIi'2a{s) H-^^^ 4 , 2 ( 5 ) - kfi{s) = 0, 

mi 


c( s') 

792 , 4 ( 5 ) + 072 ^ 2 , 4 ( 5 ) H-^^^ 4 , 2 ( 5 ) - 47 . 8 ( 5 ) = 0, 

m 2 

j^is) - 071045(5) + 7 l^^lDi^ 2 ( 5 ) = 0 , 
je{s) - uiBj&is) + 7107101 ^ 4 ( 5 ) = 0, 
j 7 (s) - 072047(5) + 72 ^^ 2 D^ 92 ( 5 ) = 0 , 
js{s) - 072048(5) + 72l^2D^94(5) = 0, 


Jl,2{s) 

-^ 3 , 4 ( 5 ) 

<^ 5 , 6 ( 5 ) 

>^ 7 , 8 ( 5 ) 

41 , 2 ( 5 ) 

43 , 4 ( 5 ) 

45 , 0 ( 5 ) 

47 , 8 ( 5 ) 


dK7i07ioe ^Vi,3(m), 

5 

du 72t^2De"‘^"°^“"^Vl,3(M), 

5 

du 7l07ioe"‘^^°^“"^V2,4(M), 
du 72072 Oe"‘^^°^“"^V 2 , 4 (M), 

5 

du 7i07ioe"‘^i°(“"^) z/ 1,3 {u), 
''*d7i72a;2De-‘^^°(“-^)79i,3(M), 

dw 7107106"^^!°1/2,4 (m) , 

5 

H 

dU72W2De““^'^'“““*?^2,4('«), 


(42) 


where the Eqs. (38-41) are solved by a Runge-Kutta algorithm. A similar process can 
be followed for arbitrary spectral densities, where the Eqs. (35-36) have to be solve by 
a nnmerical method for integro-differential eqnations. In order to give an idea, the hnal 
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equations to solve for any spectral density are 


V^i,3(s) + + 7i(0)](^i,3(s) + 

‘^2,4:{s) + [uj\ + 7l(0)](p2,4(’S) + 

01, 3 ( 5 ) + [ul + 72(O)]01,3(s) + 
02,4(’S) + [uj\ + 72(O)]02,4(’S) + 

i>i,3(s) + [ul + 7i(0)]z2i,3(s) + 
^^ 2 , 4 ( 5 ) + [ul + 7i(0)]z/2,4(’S) + 
di,3(s) + [ul + 72(0)]^9 i, 3 (-s) + 
■^2,4 (s) + + 72(0)]'d2,4(>s) + 


c(s) 

- { 

rui 

c(g) 

mi 

c{s) 

m 2 

c(g) 

m 2 

c{s) 

mi 

c{s) 

mi 

c(g) 

m 2 

c(g) 

m 2 


3,1 (-5) + j dn 7 i(s 

04,2 (-s) + / dn7i(s 
Jo 

<^ 3 ,i(-s)+ / dn72(s 

Jo 

<^ 4 , 2 ( 5 )+ [ dM72(s 


^ 3 ,i(-s) + j du-f'^iu 
■^ 4 , 2 ( 5 ) + ^ dM7^(M 
^ 3 ,i{s) + dn 72 (n 
^4,2{s) + [ dn72(n 


u)ipi,3{u) = 0 , 


u)ip2,4{u) = 0 , 


m)01,3(m) = 0, 


m)02,4(m) = 0 , 


s)z/i, 3 (n) = 0, 


s)z/ 2 , 4 (n) = 0, 


s)^l,3{u) = 0 , 


s)^2,4{u) = 0 . 


(43) 


(44) 


SQUEEZING GENERATION BY NON-MARKOVIAN DYNAMICS 

In general, for non-Markovian dynamics the system does not always thermalize and 
part of its qnantnm coherence could be preserved in the steady state’. In onr Ref. [15] 
arXiv; 1401.1418, we showed that even at thermal eqnilibrinm, non-Markovian dynamics 
allow for (i) generating sqneezing in a single harmonic mode that indnces deviations from 
the canonical thermal state and (ii) the presence of entanglement between two harmonic 
modes. For measuring cooling based on the mean phonon number this may be problematic 
because the hnal state may not be thermal. However, the presence of quantum correlations 
assisted by non-Markovian dynamics is possible at equilibrium only when k-BT/hjuj^ 1. 

For the case at hand, interest is in the the mechanical mode that is assnmed to be in a 
thermal state at t = 0 with /cbT/^ih 3> 1; therefore, there is no non-Markovianally-induced 
squeezing initially. During the subsequent ultrafast dynamics, of the order of the mechanical 
period oscillation, the time-dependent character of the conpling may sqneeze the normal 
modes of the optomechanical systems and, certainly, take the system into a non-thermal state. 
Note that this time-depedent-conpling-indnced sqneezing is also present in the Markovian 
case considered in Ref. [8] Phys. Rev. Lett. 107, 177204 (2011) and in one of the original 


10 



proposal of sideband cooling [see Refs., 7 and 22 in the mannscript]. 

A possible way to characterize the deviations from a thermal state is the g‘^{t) correlation 
fnnction; however, dne to the Ganssian character of the state of the mechanical and optical 
states, and for the present purposes, an equivalent calculation to the function is the 
direct calculation of the squeezing parameter r{t) dehned as 


r = 




(45) 


being (Ag)^ and (Ap)^ the dispersion of the position and momentum, respectively. 

For the parameters used to generate the mean phonon dynamics depicted by the black 
curve in Fig. 1 of the manuscript. Fig. f depicts the time dependence of the squeezing 
parameter of the mechanical mode state r{t). Because the initial state is thermal at high 
temperature and no initial system-bath correlations are considered, we have that r( 0 ) = 0 . 
The subsequent time-modulation of the coupling induces squeezing that goes very close 
to zero at the end of the cooling protocol. Specihcally, rNM(Gooi) = 7.f93 x 10“^ and 
?’m(Gooi) = 3.572 X 10"^. 

Therefore, non-Markovian dynamics generate more squeezing than the Markovian one; 
however, due to the parameter regime, the excess of squeezing is very small Ar(tcooi) = 
3.621 X 10“^. Thus, it is safe to consider the mean phonon number as an measure of cooling. 



FIG. 1. Squeezing parameter r{t) for the cooling protocol depicted by the black curve in Fig. 1 
of the manuscript under Markovian (blue curve) and non-Markovian (red curve) dynamics. The 
parameters are 7 = k = 2.15 x lO^^Wc, ut = 100 and Ucav = 0. 
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ENTROPY TRANSFER 


Interest concerns how to reduce the phonon number in the mechanical mode to achieve, 
as soon as possible, a state near to the ground state. An alternative cooling witness is the 
entropy at the mechanical mode. If the entropy of the mechanical mode decreases as the 
optical mode entropy increases, the effective temperature of the mechanical mode decreases, 
i.e., we are cooling the mechanical mode although we cannot give an absolute value for the 
temperature. Figure 2 depicts the von Newman entropy [2] 

+ 0 log - 0 log ~ 

X = \J {q^){p^) - {pq + qvY/^i (47) 

for the mechanical and the optical mode as a function of time. 



FIG. 2. von Neumann entropy of the both mechanical (Fm) and optical mode {Sc)- Parameters as 
in Fig. 1. Continuous lines for non-Markovian dynamics and dashed lines for Markovian dynamics. 

Cooling at this level is guaranteed by the fact that the entropy in the mechanical mode 
decreases; specihcally, it decreases from S'm(O) = 5.6101 to S'm(tcooi) = 2.23 x 10“^ in the 
Markovian case and to S'm(tcooi) = 2.19 x 10“^ in the non-Markovian case. Thus, the systems 
sits very close to its ground state (zero entropy). For the cavity mode, it varies from S'c(O) = 0 
to S'c(tcooi) = 5.601 in the Markovian case and to S'c(fcooi) = 5.607 in the non-Markovian. 
Although for this set of parameters the entropy difference is not substantial, from the analysis 
in the manuscript, it is clear that if the dissipative rate 7 increases, then the entropy difference 
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does so. Note that non-Markovian dynamics favours entropy flow as discussed in Ref. [5]. 


OPTIMIZATION ALGORITHM 

To find the minimum phonon number in the mechanical mode, Eqs. (38-41) are optimized 
to minimize the phonon number in the mechanical mode. In this algorithm the condition to 
minimize is then 

J = h[x{t),t]+ f V[x{t),g{t),t] dt = {n[x{t),t]) (48) 

Jto 

where h [x(t), f] and V [x(f), g(f), t] are arbitrary functions, x(t) and g(t) are the vector of the 
equations of motion and the vector of the optimal functions, respectively. In this work, from 
Eq. (48), R [x(f), g(f), f] = 0, h[x(f),f] = (h[x(f),f]), g(t) = c(t) being c(t) the coupling 
function and 

H = Vlx(t),g(t),t] + p-x = p-x, (49) 

where p is the vector of the costate equations. 

However, the condition to minimize, the phonon number in the mechanical mode, is given 
in terms of the second moments [see Eq. (8)] which in turn are given in terms of the functions 
ifi, (pi, Ui and 'dj. This makes quite tedious the calculation of the “initial” conditions for the 
costate equations (see below) . 

The algorithm used to solve the optimization problem was “t/ie method of steepest deseenf' 
[3] whose algorithm consists of 4 steps: 

1. Subdivide the interval [to,tf] Into N equal subintervals and assume an Initial piecewise- 
constant control g^^^t) = g{0){tk), t G [tk,tk+i] fc = 0,1, ...,N — 1. In this case, the 
control function is the coupling between the two resonators c(s) [see Eqs. (38-41)], the 
mechanical and the electromagnetic mode and the vector of the optimal functions is 
composed by one function. 

2. Apply the assumed control g^^'^ to integrate the state equations from to to tf with initial 
conditions x(to) = xq and store the state trajectory In this case, Eqs. (38-41) are 
solved and the state trajectory is stored to modify the control function later. 
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3. Apply gfW and to integrate costate eqnations backward, i.e., from [t/,to]- The 
“initial valne” is obtained by: 

rxW(t/-)l d (h f)]) 

p'-'Htf) = ^ ' ■ ( 50 ) 

Costate equations and their “initial conditions” are given below in the section “Costate 
equations and their initial conditions in the optimization algorithm (Third Step)”. 

Evaluate dH^'^\t)/dg, t G [to,tf] and store this vector. 


4. If 


dH^i) 

dg 


< e, 


(51) 


dH^i) 

^ C/ ■ 

dH^i) 

1 T r 

QH^i) 

dg 

Jto . 

dg 


dg 


(52) 


then stop the iterative procedure. Here e is a preselected small positive constant used as 
a tolerance. If Eq. (51) is not satished, adjust the piecewise-constant control function 


by: 


g0+‘>(4) = g“>-T(^(4) 

dg 


k = 0,l,...,N -1. 


(53) 


Replace g^'^'> by ^b+i) ^nd return to step 2. Here, r is the step size. 


COSTATE EQUATIONS AND THEIR INITIAL CONDITIONS IN THE 
OPTIMIZATION ALGORITHM (THIRD STEP) 


The costate equations are calculated by 


P = 


dH 
dx.' 


(54) 


The “initial” (hnal) conditions for the costate equations in Eq. (50), in terms of the 
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seconds moments, are given by 


Pviitf) = 
= 

P<Pi{if) = 

P'Pzitf) = 

Pvlpi^f) ~ 

Pip3p{tf) — 

Plpipi^f) ~ 

Pipspi^f) ~ 


1 (d{p^) 

2 V d(pi 

1 fd{p^) 

2 V d(p3 

1 (d{p^) 

2 V <901 

1 fd{p^) 

2 V <903 

1 (d{p^) 

2 \d(pip 

1 (d{p^) 

2 \d(p3p 

1 (d{p^) 

2 V 001 P 

1 (d{p^) 

2 V003P 


d{e) \ 

dpj’ 

dm ) 

9^3 J’ 

dm ) 

001 y ’ 

0(g^) \ 

003 

0(g^) \ 

0 V?iP/ ’ 

0(g^) \ 
dpspj ’ 
0(g^) \ 
00 iP/ ’ 
0(g^) \ 
003P/ ’ 


Pp2itf) = 
PpAtf) = 
P^f) = 
P^f) = 
Pv2p{tf) ~ 
P^Api^f) ~ 
P4>2pi^f) ~ 
Plp4pi^f) ~ 


If dm .dm) 

2 0 dp2 dp2 ) ’ 

1 fdm dm) 

2 0 dpi dpi ) ’ 

1 fdm dm) 

2 V 002 002 y ’ 

1 fdm dm) 

2 \ 004 004 y ’ 

1 /0(0^) 0(g^)\ 

2 V 0 <^ 2 P 0 <^ 2 Py ’ 

1 fdm dm) 

2 v0v?4P 0<^4py ’ 
ifdj^ djm) 

2 0 002P 002P y ’ 

1 ^0(0^) dm) 

2 \ 004P dcjy^p ) ’ 


(55) 


= o 


dm , 0 (g^), [dm , 0 (g^)^ . , ^ [dm , 0 (g^)^ . 

-5;7 + ^ + 1 .(a;y + j“'“ + 2 . (ay + yy)”' 


(56) 


P^f) = o 


dm , dm , ^y 0 (p^), dm\ . 

dv2 dv2 ““ \ 0^2* 0^2* y .j“ V 0z^*2 0z/i2 y 

2 — 1 2 — 1 ^ 2 y ^2 


(57) 


P^f) = o 


0 (0^), dm , ^y0(0^), dm\ fdm^dm) . 

0 Z /3 0 Z /3 \ du^i du^i ) \ duiS duiS) 

2 = 1 ^ ^ 2 = 1 ,2^3 ^ ^ 


(58) 


P^f) = o 


dm^dm ^fdm dm) ^ fdm dm) . 
0Z/4 0Z/4 I 01/4? 01/4? / I 0Z/i4 duiA J 

2 = 1 ^ ^ 2=1 2=^:4 ^ ^ 


( 59 ) 
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= o 


A (d{p^) d{f)\ 

d'&i d'di \ d'dii d-Pii) ^ .j“ \ dMl dMl) ^ 

t - 1 'h -1^2^—1 


( 60 ) 


P^2itf) = o 


d^ d^ ^(d^ d^\ A (d{p^) d{f)\ 

d'&2 d'&2 ^ V 9^2^ ) ^ ^ V d'&i‘2‘ d'&i2) ^ 

i=l ^ ^ 1 = 1 , 17^2 ^ ^ 


(61) 


P^3itf) = o 


W A (d{p^) d{f)\ 

d'ds d'&s ^ V d'dzi d'&si J I ddii d'&i?) ) 

2=1 ^ ^ 2 = 1 ,2^3 ^ ^ 


(62) 


P^A^s) = o 


^ (d^f) d{?^)\ 

dp4 d'&4 \ d'&4i d'&4i) * \ d'&iA dMA) ^ 


(63) 


with uii^r and '&ii^r given by 


uir = 
■dir = 
uii = 
'Oil = 



du Ki{u — s)ui{s), 


duK2{u - s)vi{s), 


dn Ki{u — s)ui{u), 


du K 2 {u — s)vi{u), 


(64) 


Jo Jo 

where Ki{s) is the noise kernel generated by the dissipation kernels in Eqs. (37) [1]. 
Variances 

The initial conditions for the thermal states of the mechanical (1) and optical (2) modes 
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are 


(9i)(0) = 0, 
(Pi)(0) = 0, 


(PiPi)(0) = , 

(giPi)(0) = 0, 


(g2)(0) = 0, 

(P2)(0)=0, 


( 92 ^ 2 )(0) 

(P2P2)(0) 


-coth 

27722022 


1712002 

2 


coth 



(?2P2)(0) = 0, 


where /3i are the inverse temperatures of the baths, Ui are the frequencies of the mechanical 
and electromagnetic modes respectively, and rrii the masses of the resonators. 

The second moments of the mechanical resonator {qiqi){t) = sqlql[t] and (piPi)(t) = 
splpl[t] in terms of the functions ipi = \i\, 0^ = Vi, z/* = ui and d* = vi are given by 

sqlql[t_] : = (sp2p2i (U2p [0] U4[t]-U2[t] U4p [0] ) ~2) / (4 (2 U2 [t] V2[t]-2 U4[t] 

V4[t])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t] ) “2) ~2) +(vep2i''2 (U2p[0] U4[t]-U2[t] 
U4p[0])~2)/(4 (2 U2[t] V2[t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)~2) 
+ ((U2p[0] U4[t]-U2[t] U4p[0])''2 (u22 [t] (u4[t] vl[t]-u3[t] v2 [t] ) ~2+u2 [t] “2 
(u44[t] vl [t]''2-v2 [t] (-(u33 [t]+vll [t]) v2[t]+vl[t] (u34 [t]+u43 [t]+vl2 [t]+v21 [t])) 

+ vl[t]~2 v22[t] )+u24[t] u3[t] u4[t] vl [t] v4[t]-u23[t] u4[t]''2 vl [t] v4[t]-u32[t] 
u4[t]~2 vl [t] v4[t]+u3[t] u4[t] u42 [t] vl [t] v4 [t]-u4 [t] ~2 vl [t] vl4[t] v4[t]-u24[t] 
u3[t]''2 v2 [t] v4[t]+u23[t] u3[t] u4[t] v2 [t] v4[t]+u3[t] u32[t] u4[t] v2 [t] 
v4 [t]-u3 [t] “2 u42[t] v2 [t] v4[t]+u3[t] u4[t] vl4[t] v2 [t] v4[t]+u3[t] u4[t] 
vl [t] v24[t] v4 [t]-u3 [t] ~2 v2 [t] v24[t] v4[t]-u3[t] u34[t] u4[t] v4 [t]''2+u33 [t] 
u4[t]''2 v4 [t] ~2-u3 [t] u4[t] u43[t] v4 [t]''2+u3 [t] ~2 u44[t] v4 [t] "2+u4 [t] “2 
vll[t] v4 [t] ~2-u3 [t] u4[t] vl2[t] v4 [t]''2-u3 [t] u4[t] v21 [t] v4 [t] ~2+u3 [t] "2 
v22[t] v4 [t]''2-u4 [t] "2 vl [t] v4[t] v41[t]+u3[t] u4[t] v2[t] v4[t] v41[t]+u3[t] 
u4[t] vl [t] v4[t] v42 [t]-u3 [t] ~2 v2 [t] v4[t] v42[t]+u2[t] (u24[t] vl [t] (-u4[t] 
vl[t]+u3[t] v2[t])+u4[t] (-u42[t] vl [t]''2+u23 [t] vl [t] v2[t]+u32[t] vl [t] 
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v2[t]+vl[t] vl4[t] v2 [t]-vl [t] ~2 v24 [t]+u34 [t] vl [t] v4[t]+u43[t] vl [t] v4[t]+vl[t] 
vl2[t] v4[t]-2 u33 [t] v2[t] v4[t]-2 vll[t] v2[t] v4[t]+vl[t] v21[t] v4[t]+vl[t] 
v2[t] v41 [t]-vl [t] ~2 v42 [t] )+u3 [t] (u42[t] vl [t] v2[t]-u23[t] v2 [t] ~2-u32 [t] 
v2[t]''2-vl4[t] v2[t]~2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2[t] 
v4[t]+u43[t] v2[t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] v21[t] v4[t]-2 vl [t] v22 [t] 
v4 [t]-v2 [t] “2 v41[t]+vl[t] v2 [t] v42 [t]) ) +(u4 [t] vl[t]-u3[t] v2[t])''2 v44[t]))/(2 
(u2[t] v2[t]-u4[t] v4[t])~2 (2 U2[t] V2 [t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t] )''2)''2) - (vepli vep2i (U2p[0] U4[t]-U2[t] U4p[0]) (-mU2[t] V2p[0]+m 
U4[t] V4p[0]))/(2 (2 U2[t] V2[t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)~2)-(vep2i veqli (U2p[0] U4[t]-U2[t] U4p [0]) (U2[t] (-Ulp[0] 
V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] (U2p [0] V2[t]-U4p[0] 
V4[t])) (-m U2[t] V2p[0]+m U4[t] V4p[0] ))/((2 U2[t] V2[t]-2U4[t] V4[t])~3 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 

V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t] ) “2) ~2)- (sq2p2i (U2p [0] U4[t]-U2[t] 

U4p[0]) (U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+ni U4[t] V4p[0]))/((2 U2[t] 

V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2) "2)-(vep2i veq2i 
(U2p[0] U4[t]-U2[t] U4p[0]) (U2 [t] (U4p [0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] 
Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+ni 
U4[t] V4p[0]))/((2 U2[t] V2[t]-2 U4 [t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)-2) 
- ((U2p[0] U4[t]-U2[t] U4p[0]) (ul2[t] u2 [t] u4[t] vl [t] v2[t]+u2[t] u21[t] 
u4[t] vl [t] v2 [t]-u2 [t] “2 u41 [t] vl [t] v2[t]+ul3[t] u2[t]''2 v2 [t]''2-ul2 [t] 
u2 [t] u3[t] v2 [t] ~2-u2 [t] u21[t] u3 [t] v2 [t]''2+u2 [t] "2 u31[t] v2 [t]''2+u2 [t] "2 
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vl3[t] v2 [t]''2-u2 [t] "2 vl [t] v2[t] v23[t]-2 u2 [t] u24[t] u4[t] vl [t] v3[t]+2 

u22[t] u4[t]''2 vl [t] v3[t]-2 u2[t] u4[t] u42 [t] vl [t] v3[t]+2 u2[t]''2 u44[t] 

vl [t] v3[t]+u2[t] u24[t] u3[t] v2[t] v3 [t]-u2 [t] "2 u34 [t] v2 [t] v3[t]+u2[t] 
u23[t] u4[t] v2[t] v3[t]-2 u22 [t] u3[t] u4[t] v2[t] v3[t]+u2[t] u32 [t] u4[t] 
v2[t] v3[t]+u2[t] u3[t] u42[t] v2[t] v3 [t]-u2 [t] "2 u43 [t] v2 [t] v3 [t]-u2 [t] ~2 
vl2[t] v2 [t] v3[t]+u2[t] u4[t] vl4[t] v2 [t] v3 [t]-u2 [t] ~2 v2 [t] v21[t] v3[t]+2 
u2[t]''2 vl [t] v22[t] v3[t]-2 u2[t] u4[t] vl [t] v24[t] v3[t]+u2[t] u3[t] v2[t] 
v24[t] v3 [t]+u2 [t] ~2 v2[t]"2 v31 [t]-u2 [t] ~2 vl [t] v2 [t] v32[t]+u2[t] u4[t] 
vl [t] v2[t] v34[t]-u2[t] u3[t] v2[t]''2 v34 [t]-ul2 [t] u4[t]~2 vl [t] v4[t]-u21[t] 
u4[t]~2 vl [t] v4[t]+u2[t] u4[t] u41 [t] vl [t] v4[t]-2 ul3[t] u2[t] u4[t] v2[t] 
v4[t]+ul2[t] u3[t] u4[t] v2[t] v4[t]+u21[t] u3[t] u4[t] v2[t] v4[t]-2 u2[t] 

u31[t] u4[t] v2[t] v4[t]+u2[t] u3 [t] u41 [t] v2[t] v4[t]-2 u2[t] u4 [t] vl3[t] 

v2[t] v4[t]+u2[t] u4[t] vl [t] v23[t] v4[t]+u2[t] u3 [t] v2 [t] v23[t] v4[t]+u24[t] 
u3 [t] u4[t] v3[t] v4[t]+u2[t] u34[t] u4[t] v3[t] v4[t]-u23[t] u4[t]''2 v3[t] 
v4[t]-u32[t] u4[t]~2 v3[t] v4[t]+u3[t] u4[t] u42[t] v3 [t] v4[t]+u2[t] u4[t] 
u43[t] v3[t] v4[t]-2 u2[t] u3[t] u44[t] v3[t] v4[t]+u2[t] u4[t] vl2[t] v3[t] 
v4 [t]-u4 [t] “2 vl4[t] v3[t] v4[t]+u2[t] u4[t] v21[t] v3 [t] v4[t]-2 u2 [t] u3[t] 
v22[t] v3[t] v4[t]+u3[t] u4[t] v24[t] v3[t] v4[t]-2 u2 [t] u4[t] v2 [t] v31[t] 
v4[t]+u2[t] u4[t] vl [t] v32[t] v4[t]+u2[t] u3[t] v2 [t] v32[t] v4 [t]-u4 [t] “2 
vl [t] v34[t] v4[t]+u3[t] u4[t] v2[t] v34[t] v4[t]+ul3[t] u4[t]"2 v4 [t]''2+u31 [t] 
u4[t]~2 v4 [t]''2-u3 [t] u4[t] u41 [t] v4 [t]''2+u4 [t] ~2 vl3[t] v4 [t] "2-u3 [t] u4[t] 
v23[t] v4 [t] ~2+u4 [t] "2 v31[t] v4 [t]''2-u3 [t] u4[t] v32[t] v4 [t]''2-ul4 [t] (u2[t] 
vl[t]-u3[t] v4[t]) (u2[t] v2[t]-u4[t] v4[t])+u2[t] u4[t] v2[t] v3 [t] v41 [t]-u4 [t] “2 
v3[t] v4[t] v41[t]-2 u2 [t] u4[t] vl [t] v3[t] v42[t]+u2[t] u3[t] v2 [t] v3[t] 
v42[t]+u3 [t] u4 [t] v3 [t] v4[t] v42[t]+u2 [t] u4 [t] vl [t] v2[t] v43[t]-u2 [t] 
u3 [t] v2[t]~2 v43 [t]-u4 [t] “2 vl [t] v4[t] v43[t]+u3[t] u4[t] v2[t] v4[t] v43[t]+2 
u4[t] (u4[t] vl[t]-u3[t] v2 [t]) v3[t] v44[t]+ul[t] (2 u22[t] v2 [t] (-u4[t] 
vl[t]+u3[t] v2[t])+u2[t] (u24[t] vl [t] v2[t]+u42[t] vl [t] v2[t]-u23[t] v2 [t]''2-u32 [t] 
v2[t]''2-vl4[t] v2[t]~2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2[t] 
v4[t]+u43[t] v2 [t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] v21[t] v4[t]-2 vl [t] v22 [t] 
v4[t]-v2[t] "2 v41[t]+vl[t] v2[t] v42 [t])+v4 [t] (u24[t] (u4[t] vl [t]-2 u3[t] 
v2[t])+u4[t] (-(u34[t]+u43[t] +vl2 [t]+v21 [t] ) v4[t]+v2[t] (u23 [t]+u32 [t] + 
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vl4[t] + v41[t]) + vl[t] (u42[t] + v24[t] + v42[t]))-2 u3[t] (-(u44 [t]+v22 [t]) 
v4[t]+v2[t] (u42[t] + v24[t] + v42[t])))+2 v2[t] (-u4[t] vl[t]+u3[t] v2 [t]) 
v44[t])) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 (u2 [t] v2[t]-u4[t] v4[t])~2 
(2 U2[t] V2[t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t]) ~2) ~2) + (splpli 
(-m U2[t] V2p[0]+m U4[t] V4p [0]) ~2)/(4 (2 U2[t] V2[t]-2 U4[t] V4[t])-2 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 
V2[t]-2 U4[t] V4[t] )~2)''2) + (vepli~2 (-mU2[t] V2p [0]+m U4 [t] V4p[0])''2)/(4 
(2 U2[t] V2[t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t]) ~2) ~2) + (sqlpli 
(U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0] +ni U4[t] V4p [0] ) ~2) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2)''2) + (vepli veqli 
(U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0] +ni U4[t] V4p [0] ) ~2) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2)''2) + (sqlqli 
(U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t]))~2 (-m U2 [t] V2p [0]+m U4[t] V4p [0] ) ~2) / ( (2 U2 [t] 
V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2)''2) + (veqli''2 
(U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t]))~2 (-m U2 [t] V2p [0]+m U4[t] V4p [0] ) ~2) / ( (2 U2 [t] 
V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
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(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2)''2) + (vepli veq2i 
(U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] 

(U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+ni U4[t] V4p [0] ) ~2) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2)''2) + (2 veqli 
veq2i (U2 [t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (U2[t] (U4p [0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] 
Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+ni 
U4[t] V4p[0])~2)/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p [0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)~2) + (sq2q2i (U2 [t] (U4p [0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] 
Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t]))~2 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0])~2)/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)~2) + (veq2i~2 (U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] 
Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t]))~2 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0])~2)/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 

(U2p[0] U4[t]-U2[t] U4p[0]) (V2p [0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 

U4[t] V4[t] )''2)"2) + ((-ul4[t] u2[t]''2 v2[t] v3[t]+ul2[t] u2 [t] u4[t] v2[t] 

v3[t]+u2[t] u21[t] u4[t] v2[t] v3 [t]-u2 [t] “2 u41 [t] v2 [t] v3 [t]-u2 [t] "2 v2 [t] 

v23[t] v3[t]-u2[t] u24[t] u4[t] v3 [t] "2+u22 [t] u4[t]~2 v3 [t] ~2-u2 [t] u4[t] 

u42[t] v3 [t]''2+u2 [t] “2 u44[t] v3 [t] "2+u2 [t] “2 v22 [t] v3 [t] ~2-u2 [t] u4[t] 

v24[t] v3 [t]''2-u2 [t] "2 v2[t] v3[t] v32 [t]+u2 [t] ~2 v2[t]~2 v33[t]+u2[t] u4[t] 

v2[t] v3[t] v34 [t]+ul4 [t] u2[t] u4[t] v3[t] v4[t]-ul2[t] u4[t]"2 v3[t] v4[t]-u21[t] 

u4[t]''2 v3[t] v4[t]+u2[t] u4[t] u41 [t] v3[t] v4[t]+u2[t] u4[t] v23 [t] v3[t] 

v4[t]+u2[t] u4[t] v3[t] v32[t] v4[t]-2 u2[t] u4[t] v2[t] v33[t] v4 [t]-u4 [t] ~2 

v3[t] v34[t] v4 [t]+u4 [t] ~2 v33[t] v4 [t]''2+ull [t] (u2 [t] v2[t]-u4[t] v4 [t])''2-u2 [t] 

u4[t] v3[t]~2 v42[t]+u2[t] u4[t] v2 [t] v3[t] v43 [t]-u4 [t] “2 v3[t] v4[t] v43 [t]+u4 [t] ~2 
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v3[t]''2 v44[t]+ul[t]''2 (u22[t] v2 [t]''2-v4 [t] (-(u44 [t]+v22 [t] ) v4[t]+v2[t] 

(u24 [t]+u42 [t]+v24 [t]+v42 [t] ) ) +v2[t]''2 v44 [t] )+ul [t] (ul2[t] v2[t] (-u2[t] 
v2[t]+u4[t] v4[t])+u2[t] (-u21[t] v2 [t]''2+u24 [t] v2 [t] v3[t]+u42[t] v2[t] 
v3[t]+v2[t] v24[t] v3 [t]-v2 [t] ~2 v34 [t]+ul4 [t] v2 [t] v4[t]+u41[t] v2 [t] v4[t]+v2[t] 
v23[t] v4[t]-2 u44[t] v3[t] v4[t]-2 v22[t] v3[t] v4[t]+v2[t] v32[t] v4[t]+v2[t] 
v3[t] v42 [t]-v2 [t] "2 v43 [t] )+u4 [t] (-2 u22[t] v2[t] v3[t]+v4[t] (-(ul4[t] 

+ u41[t] + v23[t] +v32[t]) v4[t]+v3[t] (u24 [t]+u42 [t]+v24 [t]+v42 [t] )+v2 [t] 

(u21 [t]+v34[t]+v43[t] ))-2 v2 [t] v3[t] v44[t]))) (-m U2 [t] V2p[0]+m U4[t] 

V4p[0])''2)/(2 (u2[t] v2[t]-u4[t] v4[t])''2 (2 U2[t] V2[t]-2 U4[t] V4[t])~2 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 

V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t] )-2) ~2) +(ni sq2p2i (U2p [0] U4[t]-U2[t] 

U4p[0])~2 (U2[t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] 
(-V2p[0] V4[t]+V2[t] V4p[0])))/((2 U2 [t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p[0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(ni (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t] )''2)~2) +(m vep2i veq2i (U2p [0] U4[t]-U2[t] U4p[0])''2 
(U2[t] (-Vlp[0] V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] 

(-V2p[0] V4[t]+V2[t] V4p[0])))/((2 U2 [t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p[0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t] )''2)~2)-(m vepli veq2i (U2p[0] U4[t]-U2[t] U4p [0]) (-m 
U2[t] V2p[0]+m U4[t] V4p[0]) (U2 [t] (-Vlp[0] V2[t]+Vl[t] V2p [0])+U4 [t] (Vlp[0] 
V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ) / ( (2 U2 [t] V2[t]-2 
U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t])-2) ~2) -(2 m veqli veq2i (U2p [0] U4[t]-U2[t] 
U4p[0]) (U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+m U4[t] V4p[0]) (U2 [t] (-Vlp[0] 
V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
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U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)-2) 
-(2 m sq2q2i (U2p [0] U4[t]-U2[t] U4p [0]) (U2[t] (U4p [0] Vl[t]-U3p[0] V2[t])+U4[t] 
(-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (-m U2 [t] 

V2p[0]+m U4[t] V4p[0]) (U2 [t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] (Vlp[0] 

V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ) / ( (2 U2 [t] V2[t]-2 
U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t]) “2) ~2) -(2 m veq2i~2 (U2p[0] U4[t]-U2[t] 
U4p[0]) (U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+m U4[t] V4p[0]) (U2 [t] (-Vlp[0] 
V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)-2) 

+ (m~2 sq2q2i (U2p [0] U4[t]-U2[t] U4p[0])~2 (U2 [t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] 
(Vlp[0] V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ~2) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])''4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2) "2) +(m''2 veq2i~2 
(U2p[0] U4[t]-U2[t] U4p[0])~2 (U2[t] (-Vlp[0] V2[t]+Vl[t] V2p [0])+U4 [t] (Vlp[0] 
V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ~2) / ( (2 U2 [t] V2[t]-2 
U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t])-2) ~2) +(m vep2i veqli (U2p [0] U4[t]-U2[t] 
U4p[0])~2 (U2[t] (V2p[0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] V4[t]+V2[t] V4p [0] )+U4 [t] 
(V3p[0] V4[t]-V3[t] V4p[0])))/((2 U2 [t] V2 [t]-2 U4[t] V4[t])~3 (((-U2p[0] 

V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t] )~2)~2)-(m sqlpli (U2p[0] U4[t]-U2[t] U4p [0]) (-m U2 [t] 

V2p[0]+m U4[t] V4p[0]) (U2 [t] (V2p [0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] 

V4[t]+V2[t] V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2 [t]-2 

U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
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U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t] )-2) ~2)- (m vepli veqli (U2p [0] U4[t]-U2[t] 
U4p[0]) (-m U2[t] V2p[0]+m U4[t] V4p[0]) (U2 [t] (V2p [0] V3[t]-V2[t] V3p [0])+U1 [t] 
(-V2p[0] V4[t]+V2[t] V4p[0])+U4[t] (V3p [0] V4[t]-V3[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2) "2) -(2 m sqlqli 
(U2p[0] U4[t]-U2[t] U4p[0]) (U2 [t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] 
V3[t]+Ulp[0] V4[t])+Ul[t] (U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+m 
U4[t] V4p[0]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0] )+Ul [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2[t]-2 U4[t] V4[t])~4 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 

V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t])-2) ~2) -(2 m veqli~2 (U2p [0] U4[t]-U2[t] 
U4p[0]) (U2[t] (-UlpEO] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+m U4[t] V4p[0]) (U2 [t] (V2p [0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)~2) 
-(2 m veqli veq2i (U2p [0] U4[t]-U2[t] U4p [0]) (U2 [t] (U4p [0] Vl[t]-U3p[0] 
V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p [0] V2[t]-U4p[0] V4[t])) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0] )+Ul [t] 

(-V2p[0] V4[t]+V2[t] V4p[0])+U4[t] (V3p [0] V4[t]-V3[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t]) "2)''2) + (2 m''2 veqli 
veq2i (U2p[0] U4[t]-U2[t] U4p[0])~2 (U2[t] (-Vlp[0] V2[t]+Vl[t] V2p [0])+U4 [t] 

(Vlp[0] V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0])) (U2 [t] (V2p [0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
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U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)-2) 
+ (m~2 sqlqli (U2p [0] U4[t]-U2[t] U4p[0])~2 (U2 [t] (V2p[0] V3[t]-V2[t] V3p [0] )+U1 [t] 
(-V2p[0] V4[t]+V2[t] V4p[0])+U4[t] (V3p [0] V4[t]-V3[t] V4p [0]) ) ~2) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])''4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2) "2) +(m''2 veqli~2 
(U2p[0] U4[t]-U2[t] U4p[0])~2 (U2[t] (V2p [0] V3[t]-V2[t] V3p [0] )+Ul [t] (-V2p[0] 
V4[t]+V2[t] V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ~2) / ( (2 U2[t] V2[t]-2 
U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t] ) ~2) ~2) ; 

splpl[t_] : = (4 sqlqli (Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] 

V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t])~2)/(2 U2[t] 

V2[t]-2 U4[t] V4[t])~2+(4 veqli''2 (Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] 

U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t])~2)/(2 
U2[t] V2[t]-2 U4[t] V4[t])''2+(8 veqli veq2i (Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] 
V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p [t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] 

V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3 [t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2+(4 
sq2q2i (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t])~2)/(2 U2[t] V2[t]-2 U4[t] V4[t])-2+(4 
veq2i''2 (U2p[t] U4[t] Vl[t]-U2[t] U4p [t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t])~2)/(2 U2[t] V2[t]-2 U4[t] V4[t])~2+(2 
(u22[t] v2[t]''2-v4[t] (-(u44[t]+v22[t] ) v4[t]+v2[t] (u24 [t]+u42 [t]+v24 [t]+v42 [t] ) ) 
+v2[t]"2 v44[t] ))/(u2[t] v2[t]-u4[t] v4[t] ) ~2+(sp2p2i (-U2p[t] U4[t]+U2[t] 

U4p[t] )~2)/(-m U2[t] V2p[0]+m U4[t] V4p[0])''2+(vep2i~2 (-U2p[t] U4[t]+U2[t] 
U4p[t])~2)/(-m U2[t] V2p[0]+m U4[t] V4p[0])~2+(2 (-U2p[t] U4[t]+U2[t] U4p[t])~2 
(u22[t] (u4[t] vl[t]-u3[t] v2 [t] )''2+u2 [t] “2 (u44[t] vl [t] ~2-v2 [t] (-(u33 [t]+vll [t] ) 
v2[t]+vl[t] (u34 [t]+u43 [t]+vl2 [t]+v21 [t] ) )+vl [t] "2 v22 [t] )+u24 [t] u3 [t] u4[t] 
vl [t] v4[t]-u23[t] u4[t]''2 vl [t] v4[t]-u32[t] u4[t]"2 vl [t] v4[t]+u3[t] u4[t] 
u42[t] vl [t] v4 [t]-u4 [t] "2 vl [t] vl4[t] v4[t]-u24[t] u3[t]~2 v2 [t] v4[t]+u23[t] 
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u3 [t] u4[t] v2[t] v4[t]+u3[t] u32 [t] u4[t] v2[t] v4 [t]-u3 [t] "2 u42 [t] v2 [t] 

v4[t]+u3[t] u4[t] vl4[t] v2[t] v4[t]+u3[t] u4[t] vl [t] v24[t] v4 [t]-u3 [t] “2 

v2[t] v24[t] v4[t]-u3[t] u34[t] u4[t] v4 [t] ~2+u33 [t] u4[t]~2 v4 [t] ~2-u3 [t] 

u4[t] u43[t] v4 [t]''2+u3 [t] “2 u44[t] v4 [t] "2+u4 [t] "2 vll[t] v4 [t] ~2-u3 [t] 

u4[t] vl2[t] v4 [t] "2-u3 [t] u4[t] v21[t] v4 [t] ~2+u3 [t] ~2 v22[t] v4 [t]''2-u4 [t] ~2 

vl [t] v4[t] v41[t]+u3[t] u4[t] v2[t] v4[t] v41[t]+u3[t] u4[t] vl [t] v4[t] 

v42 [t]-u3 [t] "2 v2[t] v4[t] v42[t]+u2[t] (u24[t] vl [t] (-u4[t] vl[t]+u3[t] 

v2[t])+u4[t] (-u42[t] vl [t] ~2+u23 [t] vl [t] v2[t]+u32[t] vl [t] v2[t]+vl[t] 

vl4[t] v2 [t]-vl [t] "2 v24 [t]+u34 [t] vl [t] v4[t]+u43[t] vl [t] v4[t]+vl[t] vl2[t] 

v4[t]-2 u33[t] v2[t] v4[t]-2 vll[t] v2 [t] v4[t]+vl[t] v21 [t] v4[t]+vl[t] 

v2[t] v41 [t]-vl [t] ~2 v42 [t] )+u3 [t] (u42[t] vl [t] v2[t]-u23[t] v2 [t] ~2-u32 [t] 

v2[t]''2-vl4[t] v2[t]''2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2[t] 

v4[t]+u43[t] v2[t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] v21[t] v4[t]-2 vl [t] v22 [t] 

v4 [t]-v2 [t] “2 v41[t]+vl[t] v2 [t] v42 [t]) ) +(u4 [t] vl[t]-u3[t] v2[t])''2 v44 [t]) ) / ( (u2 [t] 

v2[t]-u4[t] v4[t])''2 (-m U2[t] V2p [0]+m U4[t] V4p [0])''2)-(4 vep2i veqli (-U2p[t] 

U4[t]+U2[t] U4p[t]) (UlpEt] U2[t] V2[t]-Ul[t] U2p [t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] 

U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]))/((2 U2[t] V2[t]-2 U4[t] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))-(4 sq2p2i (-U2p[t] U4[t]+U2[t] U4p[t]) 

(U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p [t] V2[t]-U3p[t] 
U4[t] V4[t]+U3[t] U4p[t] V4[t]))/((2 U2[t] V2[t]-2 U4[t] V4[t]) (-m U2[t] 

V2p[0]+m U4[t] V4p[0]))-(4 vep2i veq2i (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] 

U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] 

U4[t] V4[t]+U3[t] U4p[t] V4[t]))/((2 U2[t] V2[t]-2 U4[t] V4[t]) (-m U2[t] 

V2p[0]+m U4[t] V4p[0])) + (2 (-U2p[t] U4[t]+U2[t] U4p [t]) (2 u22[t] v2 [t] (-u4[t] 
vl[t]+u3[t] v2[t])+u2[t] (u24[t] vl [t] v2[t]+u42[t] vl [t] v2[t]-u23[t] v2 [t]''2-u32 [t] 
v2[t]''2-vl4[t] v2[t]~2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2[t] 
v4[t]+u43[t] v2[t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] v21[t] v4[t]-2 vl [t] v22 [t] 
v4[t]-v2[t] "2 v41[t]+vl[t] v2[t] v42 [t])+v4 [t] (u24[t] (u4[t] vl [t]-2 u3[t] 
v2[t])+u4[t] (-(u34[t]+u43[t] +vl2 [t]+v21 [t] ) v4[t]+v2[t] (u23 [t]+u32 [t] + 
vl4[t] + v41[t]) +vl[t] (u42[t] + v24[t] + v42[t]))-2 u3[t] (-(u44 [t]+v22 [t]) 
v4[t]+v2[t] (u42[t] + v24[t] + v42[t])))+2 v2[t] (-u4[t] vl[t]+u3[t] v2 [t]) 
v44[t]))/((u2[t] v2[t]-u4[t] v4[t])''2 (-m U2[t] V2p [0]+ni U4[t] V4p[0])) + (2 
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vepli vep2i (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]))/((2 
U2[t] V2[t]-2 U4[t] V4[t])''2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(4 vep2i 
veqli (U2p[0] U4[t]-U2[t] U4p [0]) (U2p[t] V2[t]-U4p[t] V4[t]) (Ulp[t] U2[t] 
V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] 
V4[t]+Ul[t] U4p[t] V4[t]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)) +(4 sq2p2i (U2p [0] U4[t]-U2[t] U4p[0]) (U2p[t] V2[t]-U4p[t] 

V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3 [t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t]))/((2 U2 [t] V2[t]-2 U4[t] V4[t])~3 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(4 vep2i veq2i (U2p [0] U4[t]-U2[t] 
U4p[0]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] 
U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p[t] V4[t]))/((2 
U2[t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p[0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(4 vep2i 
veqli (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2 [t] (-Ulp[0] 
V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] (U2p [0] V2[t]-U4p[0] 
V4[t])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4 [t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) 
+ (4 sq2p2i (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2 [t] 

(U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p [0] 
V2[t]-U4p[0] V4[t])))/((2 U2[t] V2 [t]-2 U4[t] V4[t])-3 (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+ni U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2 [t]-2 
U4[t] V4[t])''2)) +(4 vep2i veq2i (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p [t] V2[t]-U4p[t] 
V4[t]) (U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] 
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(U2p[0] V2[t]-U4p[0] V4[t])))/((2 U2 [t] V2 [t]-2 U4[t] V4[t])~3 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(ni (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2)) -(2 (U2p [0] U4[t]-U2[t] U4p [0]) (U2p[t] V2[t]-U4p[t] 

V4[t]) (2 u22[t] v2[t] (-u4[t] vl[t]+u3[t] v2[t])+u2[t] (u24[t] vl [t] v2[t]+u42[t] 
vl [t] v2[t]-u23[t] v2 [t]''2-U32 [t] v2 [t]''2-vl4 [t] v2 [t] "2+vl [t] v2 [t] v24[t]-2 
u44[t] vl [t] v4[t]+u34[t] v2[t] v4[t]+u43[t] v2 [t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] 
v21[t] v4[t]-2 vl[t] v22[t] v4 [t]-v2 [t] "2 v41[t]+vl[t] v2 [t] v42 [t])+v4 [t] 

(u24[t] (u4[t] vl[t]-2 u3[t] v2[t])+u4[t] (-(u34 [t]+u43 [t] +vl2 [t]+v21 [t] ) 
v4[t]+v2[t] (u23 [t]+u32 [t] + vl4[t] + v41[t]) +vl [t] (u42 [t] + v24[t] + v42[t]))-2 
u3[t] (-(u44[t]+v22[t]) v4[t]+v2[t] (u42[t] + v24 [t] + v42[t])))+2 v2 [t] 

(-u4[t] vl[t]+u3[t] v2[t]) v44 [t]))/((u2 [t] v2[t]-u4[t] v4[t])~2 (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) ~2) ) + (2 (-U2p[t] 
U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (ul2[t] u2 [t] u4[t] vl [t] 
v2[t]+u2[t] u21[t] u4[t] vl [t] v2 [t]-u2 [t] “2 u41 [t] vl [t] v2[t]+ul3[t] u2[t]''2 
v2 [t]''2-ul2 [t] u2[t] u3[t] v2 [t] ~2-u2 [t] u21[t] u3[t] v2 [t]''2+u2 [t] ~2 u31[t] 
v2 [t]''2+u2 [t] ~2 vl3[t] v2 [t]''2-u2 [t] ~2 vl [t] v2[t] v23[t]-2 u2[t] u24[t] 
u4[t] vl [t] v3[t]+2 u22[t] u4[t]"2 vl [t] v3[t]-2 u2[t] u4[t] u42 [t] vl [t] 
v3[t]+2 u2[t]~2 u44[t] vl [t] v3[t]+u2[t] u24[t] u3[t] v2[t] v3 [t]-u2 [t] "2 
u34[t] v2[t] v3[t]+u2[t] u23[t] u4[t] v2[t] v3[t]-2 u22 [t] u3[t] u4[t] v2[t] 
v3[t]+u2[t] u32 [t] u4[t] v2[t] v3[t]+u2[t] u3[t] u42 [t] v2[t] v3 [t]-u2 [t] “2 
u43[t] v2[t] v3 [t]-u2 [t] "2 vl2[t] v2 [t] v3[t]+u2[t] u4[t] vl4[t] v2 [t] v3 [t]-u2 [t] “2 
v2[t] v21[t] v3[t]+2 u2[t]"2 vl [t] v22 [t] v3[t]-2 u2 [t] u4[t] vl [t] v24[t] 
v3[t]+u2[t] u3[t] v2[t] v24[t] v3 [t]+u2 [t] “2 v2[t]~2 v31 [t]-u2 [t] "2 vl [t] 
v2[t] v32[t]+u2[t] u4[t] vl [t] v2[t] v34[t]-u2[t] u3 [t] v2[t]~2 v34 [t]-ul2 [t] 
u4[t]~2 vl [t] v4[t]-u21[t] u4[t]"2 vl [t] v4[t]+u2[t] u4[t] u41 [t] vl [t] v4[t]-2 
ul3[t] u2[t] u4[t] v2 [t] v4[t]+ul2[t] u3[t] u4[t] v2 [t] v4[t]+u21[t] u3[t] 
u4[t] v2[t] v4[t]-2 u2[t] u31[t] u4[t] v2 [t] v4[t]+u2[t] u3[t] u41 [t] v2 [t] 
v4[t]-2 u2[t] u4[t] vl3[t] v2 [t] v4[t]+u2[t] u4[t] vl [t] v23[t] v4[t]+u2[t] 
u3 [t] v2[t] v23[t] v4[t]+u24[t] u3[t] u4[t] v3[t] v4[t]+u2[t] u34[t] u4[t] 



v3[t] v4[t]-u23[t] u4[t]''2 v3[t] v4[t]-u32[t] u4[t]"2 v3[t] v4[t]+u3[t] u4[t] 

u42[t] v3[t] v4[t]+u2[t] u4[t] u43[t] v3[t] v4[t]-2 u2[t] u3[t] u44[t] v3[t] 

v4[t]+u2[t] u4[t] vl2[t] v3[t] v4 [t]-u4 [t] “2 vl4[t] v3 [t] v4[t]+u2[t] u4[t] 

v21[t] v3[t] v4[t]-2 u2 [t] u3[t] v22[t] v3[t] v4[t]+u3[t] u4[t] v24[t] v3[t] 

v4[t]-2 u2[t] u4[t] v2[t] v31[t] v4[t]+u2[t] u4[t] vl [t] v32 [t] v4[t]+u2[t] 

u3 [t] v2[t] v32[t] v4 [t]-u4 [t] ~2 vl [t] v34[t] v4[t]+u3[t] u4[t] v2 [t] v34[t] 

v4[t]+ul3[t] u4[t]''2 v4 [t]''2+u31 [t] u4[t]~2 v4 [t]''2-u3 [t] u4[t] u41 [t] v4 [t]''2+u4 [t] “2 

vl3[t] v4 [t] ~2-u3 [t] u4[t] v23[t] v4 [t]''2+u4 [t] ~2 v31[t] v4 [t]''2-u3 [t] u4[t] 

v32[t] v4 [t]''2-U14 [t] (u2[t] vl[t]-u3[t] v4[t]) (u2[t] v2[t]-u4[t] v4[t])+u2[t] 

u4[t] v2[t] v3[t] v41 [t]-u4 [t]''2 v3 [t] v4[t] v41[t]-2 u2[t] u4[t] vl [t] v3[t] 

v42[t]+u2 [t] u3 [t] v2 [t] v3 [t] v42[t]+u3 [t] u4 [t] v3 [t] v4 [t] v42[t]+u2 [t] 

u4[t] vl [t] v2[t] v43[t]-u2[t] u3 [t] v2[t]~2 v43 [t]-u4 [t] “2 vl [t] v4[t] v43[t]+u3[t] 

u4[t] v2[t] v4[t] v43[t]+2 u4[t] (u4[t] vl[t]-u3[t] v2[t]) v3[t] v44[t]+ul[t] 

(2 u22 [t] v2[t] (-u4[t] vl[t]+u3[t] v2[t])+u2[t] (u24[t] vl [t] v2[t]+u42[t] 
vl [t] v2[t]-u23[t] v2 [t] ~2-u32 [t] v2 [t]''2-vl4 [t] v2 [t] "2+vl [t] v2 [t] v24[t]-2 
u44[t] vl [t] v4[t]+u34[t] v2[t] v4[t]+u43[t] v2 [t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] 
v21[t] v4[t]-2 vl[t] v22[t] v4 [t]-v2 [t] "2 v41[t]+vl[t] v2 [t] v42 [t])+v4 [t] 

(u24[t] (u4[t] vl[t]-2 u3[t] v2[t])+u4[t] (-(u34 [t]+u43 [t] +vl2 [t]+v21 [t] ) 
v4[t]+v2[t] (u23 [t]+u32 [t] + vl4[t] + v41[t]) +vl [t] (u42 [t] + v24[t] + v42[t]))-2 
u3[t] (-(u44[t]+v22[t]) v4[t]+v2[t] (u42[t] + v24 [t] + v42[t])))+2 v2 [t] 

(-u4[t] vl[t]+u3[t] v2[t]) v44[t] )))/((u2[t] v2[t]-u4[t] v4[t])~2 (2 U2 [t] 

V2[t]-2 U4[t] V4[t])''2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) -(2 sp2p2i 
(U2p[0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] 
V4[t]))/((2 U2[t] V2[t]-2 U4[t] V4[t])''2 (-mU2[t] V2p [0]+m U4[t] V4p [0] ) 

(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 

V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) -(2 vep2i~2 (U2p [0] U4[t]-U2[t] 

U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]))/((2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2 (-m U2 [t] V2p[0]+m U4[t] V4p [0]) (((-U2p[0] V2[t]+U4p[0] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
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(U2p[0] U4[t]-U2[t] U4p[0]) (V2p [0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)) -(4 (U2p[0] U4[t]-U2[t] U4p [0]) (-U2p[t] U4[t]+U2[t] U4p[t]) 

(U2p[t] V2[t]-U4p[t] V4[t]) (u22[t] (u4[t] vl[t]-u3[t] v2 [t])''2+u2 [t] “2 (u44[t] 

vl [t]''2-v2 [t] (-(u33 [t]+vll [t]) v2[t]+vl[t] (u34 [t]+u43 [t] +vl2 [t]+v21 [t]) )+vl [t] ~2 

v22[t] )+u24[t] u3[t] u4[t] vl [t] v4[t]-u23[t] u4[t]"2 vl [t] v4[t]-u32[t] 

u4[t]''2 vl [t] v4[t]+u3[t] u4[t] u42 [t] vl [t] v4 [t]-u4 [t] ~2 vl [t] vl4[t] v4[t]-u24[t] 

u3[t]''2 v2 [t] v4[t]+u23[t] u3[t] u4[t] v2 [t] v4[t]+u3[t] u32[t] u4[t] v2 [t] 

v4 [t]-u3 [t] “2 u42[t] v2 [t] v4[t]+u3[t] u4[t] vl4[t] v2 [t] v4[t]+u3[t] u4[t] 

vl [t] v24[t] v4 [t]-u3 [t] ~2 v2 [t] v24[t] v4[t]-u3[t] u34[t] u4[t] v4 [t]''2+u33 [t] 

u4[t]~2 v4 [t]''2-u3 [t] u4[t] u43[t] v4 [t]''2+u3 [t] ~2 u44[t] v4 [t] "2+u4 [t] “2 

vll[t] v4 [t]''2-u3 [t] u4[t] vl2[t] v4 [t]''2-u3 [t] u4[t] v21 [t] v4 [t] ~2+u3 [t] "2 

v22[t] v4 [t]''2-u4 [t] "2 vl [t] v4[t] v41[t]+u3[t] u4[t] v2[t] v4[t] v41[t]+u3[t] 

u4[t] vl [t] v4[t] v42 [t]-u3 [t] ~2 v2 [t] v4[t] v42[t]+u2[t] (u24[t] vl [t] (-u4[t] 

vl[t]+u3[t] v2[t])+u4[t] (-u42[t] vl [t]''2+u23 [t] vl [t] v2[t]+u32[t] vl [t] 

v2[t]+vl[t] vl4[t] v2 [t]-vl [t] ~2 v24 [t]+u34 [t] vl [t] v4[t]+u43[t] vl [t] v4[t]+vl[t] 

vl2[t] v4[t]-2 u33 [t] v2[t] v4[t]-2 vll[t] v2[t] v4[t]+vl[t] v21[t] v4[t]+vl[t] 

v2[t] v41 [t]-vl [t] "2 v42 [t] )+u3 [t] (u42[t] vl [t] v2[t]-u23[t] v2 [t] ~2-u32 [t] 

v2[t]''2-vl4[t] v2[t]~2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2[t] 

v4[t]+u43[t] v2 [t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] v21[t] v4[t]-2 vl [t] v22 [t] 

v4 [t]-v2 [t] “2 v41[t]+vl[t] v2 [t] v42 [t]) ) +(u4 [t] vl[t]-u3[t] v2[t])''2 v44 [t]) ) / ( (u2 [t] 

v2[t]-u4[t] v4[t])''2 (2 U2[t] V2[t]-2 U4[t] V4[t])~2 (-m U2[t] V2p [0]+m U4[t] 

V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) -(4 sqlpli (U2p[t] V2[t]-U4p[t] 

V4[t]) (UlpEt] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] 
V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/((2 
U2[t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) -(4 vepli 
veqli (U2p[t] V2[t]-U4p[t] V4[t]) (Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] 

U4[t] V3[t]-U2[t] U4p[t] V3[t]-UlpEt] U4Et] V4Et]+UlEt] U4pEt] V4Et]) (-m 
U2Et] V2pE0]+m U4Et] V4pE0]))/((2 U2 Et] V2Et]-2 U4Et] V4Et])~3 (((-U2pE0] 


30 



V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(ni (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 
V2[t]-2 U4[t] V4[t])~2)) -(4 vepli veq2i (U2p[t] V2[t]-U4p[t] V4[t]) (U2p[t] 

U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] 

U4[t] V4[t]+U3[t] U4p[t] V4[t]) (-m U2 [t] V2p[0]+m U4[t] V4p[0]))/((2 U2 [t] 

V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2)) -(8 sqlqli 
(U2p[t] V2[t]-U4p[t] V4[t]) (Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] 

U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (U2[t] 
(-UlpEO] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] (U2p[0] 
V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+m U4[t] V4p[0]))/((2 U2[t] V2[t]-2 
U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) -(8 veqli-2 (U2p[t] V2[t]-U4p[t] 

V4[t]) (Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] 
V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (U2[t] (-Ulp[0] V2[t]+U4p[0] 
V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] (U2p [0] V2[t]-U4p[0] V4[t])) 

(-m U2[t] V2p[0]+in U4[t] V4p[0]))/((2 U2 [t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p[0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 
V2[t]-2 U4[t] V4[t])''2)) -(8 veqli veq2i (U2p[t] V2[t]-U4p[t] V4[t]) (U2p[t] 

U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] 

U4[t] V4[t]+U3[t] U4p[t] V4[t]) (U2 [t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] 

(-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] (U2p[0] V2[t]-U4p[0] V4[t])) (-mU2[t] 

V2p[0]+m U4[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+ni U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])''2)) -(8 veqli veq2i (U2p[t] V2[t]-U4p[t] V4[t]) (Ulp[t] U2[t] 
V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] 
V4[t]+Ul[t] U4p[t] V4[t]) (U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] 
Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+m 
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U4[t] V4p[0]))/((2 U2[t] V2[t]-2 U4 [t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) 
-(8 sq2q2i (U2p[t] V2[t]-U4p[t] V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] 
U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p[t] V4[t]) (U2[t] 

(U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] 
V2[t]-U4p[0] V4[t])) (-m U2[t] V2p [0]+in U4[t] V4p[0]))/((2 U2[t] V2[t]-2 
U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(8 veq2i''2 (U2p[t] V2[t]-U4p[t] 

V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3 [t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t]) (U2[t] (U4p [0] Vl[t]-U3p[0] 
V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p [0] V2[t]-U4p[0] V4[t])) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/((2 U2 [t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p[0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])"2-(ni (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 
V2[t]-2 U4[t] V4[t])''2)) -(2 (U2p[t] V2[t]-U4p[t] V4[t]) (-2 ul [t] u22[t] 
v2[t]''2+2 u22[t] u4[t] v2[t] v3[t]+2 ul [t] u24[t] v2 [t] v4[t]-u21[t] u4[t] 
v2[t] v4[t]+2 ul [t] u42 [t] v2 [t] v4[t]+2 ul [t] v2 [t] v24[t] v4[t]-u24[t] 
u4 [t] v3 [t] v4[t]-u4 [t] u42 [t] v3 [t] v4[t]-u4 [t] v24 [t] v3 [t] v4 [t] -u4 [t] 
v2[t] v34[t] v4[t]+ul4[t] u4[t] v4 [t] ~2+u4 [t] u41 [t] v4[t]~2-2 ul [t] u44[t] 
v4[t]''2-2 ul[t] v22[t] v4[t]''2+u4[t] v23[t] v4 [t]''2+u4 [t] v32[t] v4 [t]''2+ul2 [t] 
v2[t] (u2 [t] v2[t]-u4[t] v4[t])+2 ul [t] v2[t] v4[t] v42[t]-u4[t] v3[t] v4[t] 
v42[t]-u4[t] v2 [t] v4[t] v43[t]+u2[t] (u21[t] v2 [t] "2-u24 [t] v2 [t] v3[t]-u42[t] 
v2[t] v3[t]-v2[t] v24[t] v3 [t]+v2 [t] “2 v34 [t]-ul4 [t] v2[t] v4[t]-u41[t] v2 [t] 
v4[t]-v2[t] v23[t] v4[t]+2 u44[t] v3[t] v4[t]+2 v22 [t] v3[t] v4[t]-v2[t] 
v32[t] v4[t]-v2[t] v3 [t] v42 [t]+v2 [t] "2 v43[t])+2 v2 [t] (-ul [t] v2[t]+u4[t] 
v3[t]) v44[t]) (-m U2[t] V2p [0]+m U4 [t] V4p [0] ) )/( (u2 [t] v2[t]-u4[t] v4[t])~2 
(2 U2[t] V2[t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) +(4 
m sqlpli (-U2p[t] U4[t]+U2[t] U4p[t]) (Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] 
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U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (V2p [0] 
V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])''2)) +(4 m vepli veqli (-U2p[t] U4[t]+U2[t] U4p[t]) (Ulp[t] U2[t] 
V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] 
V4[t]+Ul[t] U4p[t] V4[t]) (V2p[0] V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 
U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(4 m vepli veq2i (-U2p[t] U4[t]+U2[t] 
U4p[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t]) (V2p[0] V4[t]-V2[t] V4p[0]))/((2 
U2[t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p[0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(8 m sqlqli 
(-U2p[t] U4[t]+U2[t] U4p[t]) (UlpEt] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] 

U4[t] V3[t]-U2[t] U4p[t] V3[t]-UlpEt] U4Et] V4Et]+UlEt] U4pEt] V4Et]) (U2Et] 
C-UlpEO] V2Et]+U4pE0] V3Et])+U4Et] (-U2pE0] V3Et]+UlpE0] V4Et])+UlEt] (U2pE0] 
V2Et]-U4pE0] V4Et])) (V2pE0] V4Et]-V2Et] V4pE0]))/((2 U2Et] V2Et]-2 U4Et] 

V4Et])"4 (((-U2pE0] V2Et]+U4pE0] V4Et]) (-m U2 Et] V2pE0]+m U4Et] V4pE0]))/(2 
U2Et] V2Et]-2 U4Et] V4Et])~2-(m (U2p EO] U4Et]-U2Et] U4p EO]) (V2pE0] V4Et]-V2Et] 
V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])''2)) +(8 m veqli''2 (-U2pEt] U4Et]+U2Et] 
U4pEt]) (UlpEt] U2Et] V2Et]-UlEt] U2pEt] V2Et]+U2pEt] U4Et] V3Et]-U2Et] U4pEt] 

V3Et]-UlpEt] U4Et] V4Et]+uiEt] U4pEt] V4Et]) (U2Et] (-uipEO] V2Et]+U4pE0] 
V3Et])+U4Et] (-U2pE0] V3Et]+UlpE0] V4Et])+uiEt] (U2pE0] V2Et]-U4pE0] V4Et])) 

(V2pE0] V4Et]-V2Et] V4pE0]))/((2 U2Et] V2Et]-2 U4Et] V4Et])"4 (((-U2pE0] 
V2Et]+U4pE0] V4Et]) (-m U2 Et] V2p EO]+m U4Et] V4pE0]))/(2 U2Et] V2 Et]-2 U4Et] 
V4Et])~2-(ni (U2pE0] U4Et]-U2Et] U4p EO]) (V2pE0] V4Et]-V2Et] V4pE0]))/(2 U2 Et] 
V2Et]-2 U4Et] V4Et])~2)) +(8 m veqli veq2i (-U2pEt] U4Et]+U2Et] U4pEt]) (U2pEt] 
U4Et] VlEt]-U2Et] U4pEt] VlEt]-U2pEt] U3Et] V2Et]+U2Et] U3pEt] V2Et]-U3pEt] 

U4Et] V4Et]+U3Et] U4pEt] V4Et]) (U2 Et] (-UlpEO] V2Et]+U4pE0] V3Et])+U4Et] 

(-U2pE0] V3Et]+UlpE0] V4Et])+UlEt] (U2pE0] V2Et]-U4pE0] V4Et])) (V2p EO] V4Et]-V2Et] 
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V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4 [t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) 
+ (8 m veqli veq2i (-U2p[t] U4[t]+U2[t] U4p[t]) (Ulp[t] U2 [t] V2[t]-Ul[t] 

U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] 

U4p[t] V4[t]) (U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] 
V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p[0] V4[t]-V2[t] V4p[0]))/((2 
U2[t] V2[t]-2 U4[t] V4[t])''4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p[0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(8 m sq2q2i 
(-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] 

U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p[t] V4[t]) (U2[t] 
(U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p [0] 
V2[t]-U4p[0] V4[t])) (V2p[0] V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] 

V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(8 m veq2i''2 (-U2p[t] U4[t]+U2[t] 
U4p[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t]) (U2[t] (U4p [0] Vl[t]-U3p[0] 
V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p [0] V2[t]-U4p[0] V4[t])) 
(V2p[0] V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(ni (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 
V2[t]-2 U4[t] V4[t] )~2)) + (2 m (-U2p[t] U4[t]+U2[t] U4p [t]) (-2 ul [t] u22[t] 
v2[t]''2+2 u22[t] u4[t] v2[t] v3[t]+2 ul [t] u24[t] v2 [t] v4[t]-u21[t] u4[t] 
v2[t] v4[t]+2 ul[t] u42[t] v2 [t] v4[t]+2 ul [t] v2 [t] v24[t] v4[t]-u24[t] 
u4 [t] v3 [t] v4 [t] -u4 [t] u42 [t] v3 [t] v4 [t] -u4 [t] v24 [t] v3 [t] v4 [t] -u4 [t] 
v2[t] v34[t] v4[t]+ul4[t] u4[t] v4 [t] ~2+u4 [t] u41 [t] v4[t]~2-2 ul [t] u44[t] 
v4[t]''2-2 ul [t] v22[t] v4 [t]''2+u4 [t] v23[t] v4 [t]''2+u4 [t] v32[t] v4 [t]''2+ul2 [t] 
v2[t] (u2 [t] v2[t]-u4[t] v4[t])+2 ul [t] v2[t] v4[t] v42[t]-u4[t] v3[t] v4[t] 
v42[t]-u4[t] v2 [t] v4[t] v43[t]+u2[t] (u21[t] v2 [t]''2-u24 [t] v2 [t] v3[t]-u42[t] 
v2[t] v3[t]-v2[t] v24[t] v3 [t]+v2 [t] “2 v34 [t]-ul4 [t] v2[t] v4[t]-u41[t] v2 [t] 
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v4[t]-v2[t] v23[t] v4[t]+2 u44[t] v3[t] v4[t]+2 v22 [t] v3[t] v4[t]-v2[t] 
v32[t] v4[t]-v2[t] v3 [t] v42 [t]+v2 [t] "2 v43[t])+2 v2 [t] (-ul [t] v2[t]+u4[t] 
v3[t]) v44[t]) (V2p[0] V4[t]-V2[t] V4p [0] ) ) / ( (u2 [t] v2[t]-u4[t] v4[t])''2 
(2 U2[t] V2[t]-2 U4[t] V4[t])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t]) ~2)) +(2 
m sp2p2i (U2p[0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t])~2 (V2p[0] 
V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~2 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t] )-2) ) + (2 m vep2i''2 (U2p [0] U4[t]-U2[t] 

U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t])-2 (V2p [0] V4[t]-V2[t] V4p[0]))/((2 U2[t] 
V2[t]-2 U4[t] V4[t])~2 (-m U2 [t] V2p [0]+m U4[t] V4p[0])~2 (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+ni U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)) +(4 m (U2p [0] U4[t]-U2[t] U4p [0]) (-U2p[t] U4[t]+U2[t] U4p[t])~2 
(u22[t] (u4[t] vl[t]-u3[t] v2 [t] )''2+u2 [t] “2 (u44[t] vl [t] ~2-v2 [t] (-(u33 [t]+vll [t] ) 
v2[t]+vl[t] (u34 [t]+u43 [t] +vl2 [t]+v21 [t] ) )+vl [t] ~2 v22 [t] )+u24 [t] u3 [t] 
u4[t] vl [t] v4[t]-u23[t] u4[t]~2 vl [t] v4[t]-u32[t] u4[t]''2 vl [t] v4[t]+u3[t] 
u4[t] u42[t] vl [t] v4 [t]-u4 [t]''2 vl [t] vl4[t] v4[t]-u24[t] u3[t]''2 v2 [t] 
v4[t]+u23[t] u3[t] u4[t] v2[t] v4[t]+u3[t] u32 [t] u4[t] v2[t] v4 [t]-u3 [t] “2 
u42[t] v2[t] v4[t]+u3[t] u4[t] vl4[t] v2[t] v4[t]+u3[t] u4[t] vl [t] v24[t] 
v4 [t]-u3 [t] “2 v2[t] v24[t] v4[t]-u3[t] u34[t] u4[t] v4 [t] "2+u33 [t] u4[t]~2 
v4 [t]''2-u3 [t] u4[t] u43[t] v4 [t]''2+u3 [t] "2 u44[t] v4 [t] ~2+u4 [t] "2 vll[t] 
v4 [t]''2-u3 [t] u4[t] vl2[t] v4 [t]''2-u3 [t] u4[t] v21[t] v4 [t]''2+u3 [t] ~2 v22[t] 
v4 [t]''2-u4 [t] ~2 vl [t] v4[t] v41[t]+u3[t] u4[t] v2 [t] v4[t] v41[t]+u3[t] u4[t] 
vl [t] v4[t] v42 [t]-u3 [t] ~2 v2 [t] v4[t] v42[t]+u2[t] (u24[t] vl [t] (-u4[t] 
vl[t]+u3[t] v2[t])+u4[t] (-u42[t] vl [t]''2+u23 [t] vl [t] v2[t]+u32[t] vl [t] 
v2[t]+vl[t] vl4[t] v2 [t]-vl [t] ~2 v24 [t]+u34 [t] vl [t] v4[t]+u43[t] vl [t] v4[t]+vl[t] 
vl2[t] v4[t]-2 u33 [t] v2[t] v4[t]-2 vll[t] v2[t] v4[t]+vl[t] v21[t] v4[t]+vl[t] 
v2[t] v41 [t]-vl [t] ~2 v42 [t] )+u3 [t] (u42[t] vl [t] v2[t]-u23[t] v2 [t] ~2-u32 [t] 
v2[t]''2-vl4[t] v2[t]''2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2[t] 
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v4[t]+u43[t] v2[t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] v21[t] v4[t]-2 vl [t] v22 [t] 
v4[t]-v2[t] "2 v41[t]+vl[t] v2[t] v42 [t]) ) +(u4 [t] vl[t]-u3[t] v2[t])''2 v44[t]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/((u2[t] v2[t]-u4[t] v4[t])''2 (2 U2[t] V2[t]-2 
U4[t] V4[t])~2 (-m U2[t] V2p [0]+m U4 [t] V4p[0])~2 (((-U2p[0] V2[t]+U4p[0] 

V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])''2)) -(2 m vepli vep2i (-U2p[t] U4[t]+U2[t] U4p[t])''2 (V2p[0] 
V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~2 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(4 m vep2i veqli (U2p [0] U4[t]-U2[t] 
U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (Ulp[t] U2 [t] V2[t]-Ul[t] U2p [t] V2[t]+U2p[t] 
U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (V2p [0] 
V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~3 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) -(4 m sq2p2i (U2p [0] U4[t]-U2[t] 

U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p [t] Vl[t]-U2p[t] 
U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p[t] V4[t]) (V2p [0] 
V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])~3 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(4 m vep2i veq2i (U2p [0] U4[t]-U2[t] 
U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p [t] Vl[t]-U2p[t] 
U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p[t] V4[t]) (V2p [0] 
V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] V4[t])''3 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(4 m vep2i veqli (-U2p[t] U4[t]+U2[t] 
U4p[t])~2 (U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] 
V4[t])+Ul[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p[0] V4[t]-V2[t] V4p[0]))/((2 
U2[t] V2[t]-2 U4[t] V4[t])~3 (-m U2 [t] V2p [0]+m U4 [t] V4p [0]) (((-U2p[0] 
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V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2-(ni (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 
V2[t]-2 U4[t] V4[t])~2)) -(4 m sq2p2i (-U2p[t] U4[t]+U2[t] U4p[t])''2 (U2[t] 

(U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p [0] 
V2[t]-U4p[0] V4[t])) (V2p[0] V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2 U4[t] 

V4[t])~3 (-m U2[t] V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4 [t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) 
-(4 m vep2i veq2i (-U2p[t] U4[t]+U2[t] U4p[t])''2 (U2[t] (U4p [0] Vl[t]-U3p[0] 
V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p [0] V2[t]-U4p[0] V4[t])) 
(V2p[0] V4[t]-V2[t] V4p[0]))/((2 U2[t] V2[t]-2U4[t] V4[t])''3 (-mU2[t] V2p [0]+m 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4 [t] ) “2) ) + (2 m (U2p [0] U4[t]-U2[t] U4p [0] ) 
(-U2p[t] U4[t]+U2[t] U4p[t]) (2 u22[t] v2[t] (-u4[t] vl[t]+u3[t] v2[t])+u2[t] 
(u24[t] vl [t] v2[t]+u42[t] vl [t] v2[t]-u23[t] v2 [t]''2-u32 [t] v2 [t] ~2-vl4 [t] 
v2[t]''2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2 [t] v4[t]+u43[t] 
v2[t] v4[t]+vl2[t] v2[t] v4[t]+v2[t] v21 [t] v4[t]-2 vl [t] v22[t] v4 [t]-v2 [t] "2 
v41[t]+vl[t] v2 [t] v42 [t] )+v4 [t] (u24 [t] (u4[t] vl[t]-2 u3[t] v2[t])+u4[t] 

(-(u34 [t]+u43 [t] +vl2 [t]+v21 [t] ) v4[t]+v2[t] (u23 [t]+u32 [t] + vl4[t] + v41[t]) 
+vl[t] (u42[t] + v24[t] + v42[t]))-2 u3[t] (-(u44 [t]+v22 [t] ) v4[t]+v2[t] 

(u42[t] + v24[t] + v42[t])))+2 v2[t] (-u4[t] vl[t]+u3[t] v2 [t]) v44[t]) (V2p [0] 
V4[t]-V2[t] V4p[0]))/((u2[t] v2[t]-u4[t] v4[t])~2 (2 U2[t] V2[t]-2 U4[t] 

V4[t])~2 (-m U2[t] V2p[0]+m U4[t] V4p [0]) (((-U2p[0] V2[t]+U4p[0] V4 [t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) 
-(2 m (-U2p[t] U4[t]+U2[t] U4p[t])~2 (ul2[t] u2 [t] u4[t] vl [t] v2[t]+u2[t] 
u21[t] u4[t] vl [t] v2 [t]-u2 [t] ~2 u41 [t] vl [t] v2[t]+ul3[t] u2[t]~2 v2 [t]''2-ul2 [t] 
u2 [t] u3[t] v2 [t] ~2-u2 [t] u21[t] u3 [t] v2 [t]''2+u2 [t] "2 u31[t] v2 [t]''2+u2 [t] "2 
vl3[t] v2 [t]''2-u2 [t] "2 vl [t] v2[t] v23[t]-2 u2 [t] u24[t] u4[t] vl [t] v3[t]+2 
u22[t] u4[t]''2 vl [t] v3[t]-2 u2[t] u4[t] u42 [t] vl [t] v3[t]+2 u2[t]~2 u44[t] 
vl [t] v3[t]+u2[t] u24[t] u3[t] v2[t] v3 [t]-u2 [t] "2 u34 [t] v2 [t] v3[t]+u2[t] 
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u23[t] u4[t] v2[t] v3[t]-2 u22 [t] u3[t] u4[t] v2[t] v3[t]+u2[t] u32 [t] u4[t] 
v2[t] v3[t]+u2[t] u3[t] u42[t] v2[t] v3 [t]-u2 [t] "2 u43 [t] v2 [t] v3 [t]-u2 [t] ~2 
vl2[t] v2 [t] v3[t]+u2[t] u4[t] vl4[t] v2 [t] v3 [t]-u2 [t] ~2 v2 [t] v21[t] v3[t]+2 
u2[t]''2 vl [t] v22[t] v3[t]-2 u2[t] u4[t] vl [t] v24[t] v3[t]+u2[t] u3[t] v2[t] 
v24[t] v3 [t]+u2 [t] ~2 v2[t]"2 v31 [t]-u2 [t] ~2 vl [t] v2 [t] v32[t]+u2[t] u4[t] 
vl [t] v2[t] v34[t]-u2[t] u3[t] v2[t]''2 v34 [t]-ul2 [t] u4[t]~2 vl [t] v4[t]-u21[t] 
u4[t]''2 vl [t] v4[t]+u2[t] u4[t] u41 [t] vl [t] v4[t]-2 ul3[t] u2[t] u4[t] v2[t] 
v4[t]+ul2[t] u3[t] u4[t] v2[t] v4[t]+u21[t] u3[t] u4[t] v2[t] v4[t]-2 u2[t] 
u31[t] u4[t] v2[t] v4[t]+u2[t] u3 [t] u41 [t] v2[t] v4[t]-2 u2[t] u4 [t] vl3[t] 
v2[t] v4[t]+u2[t] u4[t] vl [t] v23[t] v4[t]+u2[t] u3 [t] v2 [t] v23[t] v4[t]+u24[t] 
u3 [t] u4[t] v3[t] v4[t]+u2[t] u34[t] u4[t] v3[t] v4[t]-u23[t] u4[t]''2 v3[t] 
v4[t]-u32[t] u4[t]~2 v3[t] v4[t]+u3[t] u4[t] u42[t] v3 [t] v4[t]+u2[t] u4[t] 
u43[t] v3[t] v4[t]-2 u2[t] u3[t] u44[t] v3[t] v4[t]+u2[t] u4[t] vl2[t] v3[t] 
v4 [t]-u4 [t] “2 vl4[t] v3[t] v4[t]+u2[t] u4[t] v21[t] v3 [t] v4[t]-2 u2 [t] u3[t] 
v22[t] v3[t] v4[t]+u3[t] u4[t] v24[t] v3[t] v4[t]-2 u2 [t] u4[t] v2 [t] v31[t] 
v4[t]+u2[t] u4[t] vl [t] v32[t] v4[t]+u2[t] u3[t] v2[t] v32[t] v4 [t]-u4 [t] “2 
vl [t] v34[t] v4[t]+u3[t] u4[t] v2[t] v34[t] v4[t]+ul3[t] u4[t]"2 v4 [t]''2+u31 [t] 
u4[t]~2 v4 [t] ~2-u3 [t] u4[t] u41 [t] v4 [t]''2+u4 [t] ~2 vl3[t] v4 [t] "2-u3 [t] u4[t] 
v23[t] v4 [t] ~2+u4 [t] "2 v31[t] v4 [t]''2-u3 [t] u4[t] v32[t] v4 [t]''2-ul4 [t] (u2[t] 
vl[t]-u3[t] v4[t]) (u2[t] v2[t]-u4[t] v4[t])+u2[t] u4[t] v2[t] v3 [t] v41 [t]-u4 [t] “2 
v3[t] v4[t] v41[t]-2 u2 [t] u4[t] vl [t] v3[t] v42[t]+u2[t] u3[t] v2 [t] v3[t] 
v42[t]+u3 [t] u4 [t] v3 [t] v4[t] v42[t]+u2 [t] u4 [t] vl [t] v2[t] v43[t]-u2 [t] 
u3 [t] v2[t]~2 v43 [t]-u4 [t] “2 vl [t] v4[t] v43[t]+u3[t] u4[t] v2[t] v4[t] v43[t]+2 
u4[t] (u4[t] vl[t]-u3[t] v2 [t]) v3[t] v44[t]+ul[t] (2 u22[t] v2 [t] (-u4[t] 
vl[t]+u3[t] v2[t])+u2[t] (u24[t] vl [t] v2[t]+u42[t] vl [t] v2[t]-u23[t] v2 [t]''2-u32 [t] 
v2[t]''2-vl4[t] v2[t]~2+vl[t] v2[t] v24[t]-2 u44[t] vl [t] v4[t]+u34[t] v2[t] 
v4[t]+u43[t] v2[t] v4[t]+vl2[t] v2 [t] v4[t]+v2[t] v21[t] v4[t]-2 vl [t] v22 [t] 
v4[t]-v2[t] "2 v41[t]+vl[t] v2[t] v42 [t])+v4 [t] (u24[t] (u4[t] vl [t]-2 u3[t] 
v2[t])+u4[t] (-(u34[t]+u43[t] +vl2 [t]+v21 [t] ) v4[t]+v2[t] (u23 [t]+u32 [t] + 
vl4[t] + v41[t]) +vl[t] (u42[t] + v24[t] + v42[t]))-2 u3[t] (-(u44 [t]+v22 [t]) 
v4[t]+v2[t] (u42[t] + v24[t] + v42[t])))+2 v2[t] (-u4[t] vl[t]+u3[t] v2 [t]) 
v44[t])) (V2p[0] V4[t]-V2[t] V4p [0] ) ) / ( (u2 [t] v2[t]-u4[t] v4[t])~2 (2 U2[t] 
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V2[t]-2 U4[t] V4[t])~2 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+ni U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)) +(4 sqlqlEt] (2 U2 [t] V2[t]-2 U4 [t] V4[t])~2 (-(((U2p[t] 
V2[t]-U4p[t] V4[t]) (-m U2 [t] V2p[0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])~2) + (m (-U2p[t] U4[t]+U2[t] U4p [t]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])''2)''2) /(-m U2[t] V2p [0]+m U4[t] V4p[0])"2+(4 m 
sq2p2i (-U2p[t] U4[t]+U2[t] U4p[t])~2 (U2 [t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] 
(Vlp[0] V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t]) (-m U2[t] V2p [0]+m U4[t] V4p [0])''2) + (4 m vep2i veq2i 
(-U2p[t] U4[t]+U2[t] U4p[t])~2 (U2 [t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] 

(Vlp[0] V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t]) (-m U2[t] V2p [0]+m U4[t] V4p [0] ) "2) - (8 m veqli veq2i 
(-U2p[t] U4[t]+U2[t] U4p[t]) (UlpEt] U2Et] V2Et]-UlEt] U2pEt] V2Et]+U2pEt] 

U4Et] V3Et]-U2Et] U4pEt] V3Et]-UlpEt] U4Et] V4Et]+UlEt] U4pEt] V4Et]) (U2Et] 
C-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0] )+U3Et] (-V2pE0] 
V4Et]+V2Et] V4pE0])))/((2 U2Et] V2 Et]-2 U4Et] V4Et])''2 (-m U2Et] V2p EO]+ni 
U4Et] V4pE0]))-(8 m sq2q2i (-U2pEt] U4Et]+U2Et] U4pEt]) (U2pEt] U4 Et] VlEt]-U2Et] 
U4pEt] VlEt]-U2pEt] U3Et] V2Et]+U2Et] U3pEt] V2Et]-U3pEt] U4Et] V4Et]+U3Et] 

U4pEt] V4Et]) (U2Et] C-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] 
V4pE0])+U3Et] (-V2pE0] V4Et]+V2Et] V4p EO] ) ) ) / ( (2 U2 Et] V2 Et]-2 U4Et] V4Et])~2 
(-m U2Et] V2pE0]+m U4Et] V4pE0]))-(8 m veq2i''2 (-U2pEt] U4Et]+U2Et] U4pEt]) 

(U2pEt] U4Et] VlEt]-U2Et] U4pEt] VlEt]-U2pEt] U3Et] V2Et]+U2Et] U3pEt] V2Et]-U3pEt] 
U4Et] V4Et]+U3Et] U4pEt] V4Et]) (U2Et] (-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] 

(VlpEO] V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] V4Et]+V2Et] V4pE0])))/((2 U2Et] 

V2Et]-2 U4Et] V4Et])~2 (-m U2 Et] V2pE0]+m U4Et] V4pE0])) + (4 m vepli veq2i 
(-U2pEt] U4Et]+U2Et] U4pEt]) (U2pEt] V2Et]-U4pEt] V4Et]) (U2Et] (-VlpEO] 

V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] V4Et]+V2Et] 
V4pE0])))/((2 U2Et] V2Et]-2 U4Et] V4Et])~3 (((-U2pE0] V2Et]+U4pE0] V4Et]) 

(-m U2Et] V2pE0]+in U4Et] V4pE0]))/(2 U2 Et] V2Et]-2 U4Et] V4Et])~2-(m (U2p EO] 
U4Et]-U2Et] U4pE0]) (V2pE0] V4Et]-V2Et] V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])-2)) 
+ (8 m veqli veq2i (U2p EO] U4Et]-U2Et] U4p EO]) (U2pEt] V2Et]-U4pEt] V4Et]) 
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(Ulp[t] U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p [t] V3[t]-Ulp[t] 
U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (U2[t] (-Vlp[0] V2[t]+Vl[t] V2p [0])+U4 [t] 

(Vlp[0] V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 

(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(8 m sq2q2i 
(U2p[0] U4[t]-U2[t] U4p[0]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] 
U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] 

U4p[t] V4[t]) (U2[t] (-VlpEO] V2[t]+Vl[t] V2p [0])+U4 [t] (Vlp[0] V4[t]-Vl[t] 
V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0] ) ) ) / ( (2 U2 [t] V2 [t]-2 U4[t] V4[t])~4 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 

V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(8 m veq2i~2 (U2p [0] U4[t]-U2[t] 

U4p[0]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] 
U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p[t] V4[t]) (U2[t] 

C-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0] )+U3Et] (-V2pE0] 
V4Et]+V2Et] V4pE0])))/((2 U2Et] V2 Et]-2 U4Et] V4Et])~4 (((-U2pE0] V2Et]+U4pE0] 

V4Et]) (-m U2Et] V2pE0]+m U4Et] V4pE0]))/(2 U2 Et] V2Et]-2 U4Et] V4 Et] ) "2-(m 
(U2pE0] U4Et]-U2Et] U4pE0]) (V2pE0] V4Et]-V2Et] V4pE0]))/(2 U2 Et] V2Et]-2 
U4Et] V4Et])''2)) +(8 m veqli veq2i (-U2pEt] U4Et]+U2Et] U4pEt]) (U2pEt] V2Et]-U4pEt] 
V4Et]) (U2Et] C-UlpEO] V2Et]+U4pE0] V3Et])+U4Et] (-U2pE0] V3Et]+UlpE0] V4Et])+UlEt] 
(U2pE0] V2Et]-U4pE0] V4Et])) (U2Et] (-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] 
V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] V4Et]+V2Et] V4pE0])))/((2 U2Et] V2Et]-2 
U4Et] V4Et])~4 (((-U2pE0] V2Et]+U4pE0] V4Et]) (-m U2 Et] V2p EO]+m U4Et] V4pE0]))/(2 
U2Et] V2Et]-2 U4Et] V4Et])~2-(m (U2p EO] U4Et]-U2Et] U4p EO]) (V2pE0] V4Et]-V2Et] 
V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])''2)) +(8 m sq2q2i (-U2pEt] U4Et]+U2Et] 

U4pEt]) (U2pEt] V2Et]-U4pEt] V4Et]) (U2Et] (U4pE0] VlEt]-U3pE0] V2Et])+U4Et] 

(-U2pE0] VlEt]+U3pE0] V4Et])+U3Et] (U2pE0] V2Et]-U4pE0] V4Et])) (U2 Et] (-VlpEO] 
V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] V4Et]+V2Et] 
V4pE0])))/((2 U2Et] V2Et]-2 U4Et] V4Et])~4 (((-U2pE0] V2Et]+U4pE0] V4Et]) 

(-m U2Et] V2pE0]+m U4Et] V4pE0]))/(2 U2 Et] V2Et]-2 U4Et] V4 Et] ) ~2-(m (U2p EO] 
U4Et]-U2Et] U4pE0]) (V2pE0] V4Et]-V2Et] V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])~2)) 
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+ (8 m veq2i~2 (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2 [t] 

(U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] 
V2[t]-U4p[0] V4[t])) (U2[t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] (Vlp[0] V4[t]-Vl[t] 
V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0] ) ) ) / ( (2 U2 [t] V2[t]-2 U4[t] V4[t])~4 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) -(3 m sq2p2i (U2p [0] U4[t]-U2[t] 

U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2[t] (-Vlp[0] 
V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])"3 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(3 m vep2i veq2i (U2p [0] U4[t]-U2[t] 
U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2[t] (-Vlp[0] 
V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])"3 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(3 m''2 sq2p2i (U2p [0] U4[t]-U2[t] 
U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t])-2 (V2p [0] V4[t]-V2[t] V4p[0]) (U2[t] 

(-VlpEO] V2[t]+Vl[t] V2p[0])+U4[t] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] 
V4Et]+V2Et] V4pE0])))/((2 U2Et] V2 Et]-2 U4Et] V4Et])''3 (-m U2Et] V2p EO]+m 
U4Et] V4pE0])~2 (((-U2pE0] V2Et]+U4pE0] V4Et]) (-mU2Et] V2p EO]+m U4 Et] V4pE0]))/(2 
U2Et] V2Et]-2 U4Et] V4Et])~2-(m (U2p EO] U4Et]-U2Et] U4p EO]) (V2pE0] V4Et]-V2Et] 
V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])''2)) +(3 m''2 vep2i veq2i (U2p EO] U4Et]-U2Et] 
U4pE0]) (-U2pEt] U4Et]+U2Et] U4pEt])-2 (V2p EO] V4Et]-V2Et] V4p EO]) (U2Et] 

(-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] 
V4Et]+V2Et] V4pE0])))/((2 U2Et] V2 Et]-2 U4Et] V4Et])''3 (-m U2Et] V2p EO]+m 
U4Et] V4pE0])~2 (((-U2pE0] V2Et]+U4pE0] V4Et]) (-mU2Et] V2p EO]+m U4 Et] V4pE0]))/(2 
U2Et] V2Et]-2 U4Et] V4Et])~2-(m (U2p EO] U4Et]-U2Et] U4p EO]) (V2pE0] V4Et]-V2Et] 
V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])''2)) -(4 m''2 vepli veq2i (-U2pEt] U4Et]+U2Et] 
U4pEt])~2 (V2pE0] V4Et]-V2Et] V4p EO]) (U2 Et] (-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] 
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(Vlp[0] V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~3 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])''2)) -(8 m~2 veqli veq2i (U2p [0] U4[t]-U2[t] U4p [0]) (-U2p[t] 
U4[t]+U2[t] U4p[t]) (UlpEt] U2[t] V2[t]-Ul[t] U2p [t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] 
U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (V2p[0] V4[t]-V2[t] V4p [0]) 
(U2[t] (-Vlp[0] V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] 

(-V2p[0] V4[t]+V2[t] V4p[0] )))/((2 U2 [t] V2[t]-2 U4[t] V4[t])''4 (-mU2[t] 

V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) -(8 m~2 sq2q2i 
(U2p[0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] U4[t] Vl[t]-U2[t] 
U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] 

U4p[t] V4[t]) (V2p[0] V4[t]-V2[t] V4p [0]) (U2[t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] 
(Vlp[0] V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~4 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)) -(8 m~2 veq2i~2 (U2p [0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] 
U4p[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p[t] 
V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p [t] V4[t]) (V2p[0] V4[t]-V2[t] V4p [0]) 

(U2[t] (-Vlp[0] V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] 

(-V2p[0] V4[t]+V2[t] V4p[0] )))/((2 U2 [t] V2[t]-2 U4[t] V4[t])''4 (-mU2[t] 

V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) -(8 m''2 veqli 
veq2i (-U2p[t] U4[t]+U2[t] U4p[t])"2 (U2 [t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] 
(-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p[0] V4[t]-V2[t] 
V4p[0]) (U2[t] C-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] 
(-V2pE0] V4Et]+V2Et] V4pE0] )))/((2 U2 Et] V2Et]-2 U4Et] V4Et])''4 (-mU2Et] 

V2pE0]+m U4Et] V4pE0]) (((-U2pE0] V2Et]+U4pE0] V4Et]) (-m U2 Et] V2p EO]+m 
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U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) -(8 m~2 sq2q2i 
(-U2p[t] U4[t]+U2[t] U4p[t])~2 (U2 [t] (U4p [0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] 
Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p [0] V4[t]-V2[t] 

V4p[0]) (U2[t] (-VlpEO] V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0])+U3 [t] 
(-V2p[0] V4[t]+V2[t] V4p[0] )))/((2 U2 [t] V2[t]-2 U4[t] V4[t])''4 (-mU2[t] 

V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) -(8 m''2 veq2i~2 
(-U2p[t] U4[t]+U2[t] U4p[t])~2 (U2 [t] (U4p [0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] 
Vl[t]+U3p[0] V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p [0] V4[t]-V2[t] 

V4p[0]) (U2[t] C-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] 
(-V2pE0] V4Et]+V2Et] V4pE0] )))/((2 U2 Et] V2Et]-2 U4Et] V4Et])''4 (-mU2Et] 

V2pE0]+m U4Et] V4pE0]) (((-U2pE0] V2Et]+U4pE0] V4Et]) (-m U2 Et] V2p EO]+m 
U4Et] V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])-2-(m (U2p EO] U4Et]-U2Et] U4pE0]) 
(V2pE0] V4Et]-V2Et] V4pE0]))/(2 U2 Et] V2Et]-2 U4Et] V4Et])''2)) +(4 m''2 sq2q2i 
(-U2pEt] U4Et]+U2Et] U4pEt])~2 (U2 Et] (-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] 

(VlpEO] V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] V4Et]+V2Et] V4pE0]))~2)/((2 U2Et] 

V2Et]-2 U4Et] V4Et])~2 (-m U2 Et] V2p EO]+m U4Et] V4p EO] ) ~2) + (4 m''2 veq2i~2 
(-U2pEt] U4Et]+U2Et] U4pEt])~2 (U2 Et] (-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] 

(VlpEO] V4Et]-VlEt] V4pE0])+U3Et] (-V2pE0] V4Et]+V2Et] V4pE0]))~2)/((2 U2Et] 

V2Et]-2 U4Et] V4Et])''2 (-m U2 Et] V2p EO]+m U4Et] V4p EO] ) ~2) - (8 m''2 sq2q2i 
(U2pE0] U4Et]-U2Et] U4pE0]) (-U2pEt] U4Et]+U2Et] U4pEt]) (U2pEt] V2Et]-U4pEt] 

V4Et]) (U2Et] (-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] 
(-V2pE0] V4Et]+V2Et] V4pE0]))~2)/((2 U2Et] V2Et]-2U4Et] V4Et])~4 (-mU2Et] 

V2pE0]+m U4Et] V4pE0]) (((-U2pE0] V2Et]+U4pE0] V4Et]) (-m U2 Et] V2p EO]+m 
U4Et] V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])~2-(m (U2p EO] U4Et]-U2Et] U4pE0]) 
(V2pE0] V4Et]-V2Et] V4pE0]))/(2 U2 Et] V2Et]-2 U4Et] V4Et])''2)) -(8 m''2 veq2i''2 
(U2pE0] U4Et]-U2Et] U4pE0]) (-U2pEt] U4Et]+U2Et] U4pEt]) (U2pEt] V2Et]-U4pEt] 

V4Et]) (U2Et] (-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] 
(-V2pE0] V4Et]+V2Et] V4pE0]))~2)/((2 U2Et] V2Et]-2U4Et] V4Et])~4 (-mU2Et] 

V2pE0]+m U4Et] V4pE0]) (((-U2pE0] V2Et]+U4pE0] V4Et]) (-m U2 Et] V2p EO]+m 
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U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(8 ffi-S sq2q2i 
(U2p[0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t])~2 (V2p[0] V4[t]-V2[t] 
V4p[0]) (U2[t] (-VlpEO] V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0])+U3 [t] 
(-V2p[0] V4[t]+V2[t] V4p[0]))~2)/((2 U2[t] V2[t]-2U4[t] V4[t])~4 (-mU2[t] 

V2p[0]+m U4[t] V4p[0])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2[t] V2p[0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(8 ffi-S veq2i''2 
(U2p[0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t])~2 (V2p[0] V4[t]-V2[t] 
V4p[0]) (U2[t] C-VlpEO] V2Et]+VlEt] V2pE0])+U4Et] (VlpEO] V4Et]-VlEt] V4pE0])+U3Et] 
(-V2pE0] V4Et]+V2Et] V4pE0]))~2)/((2 U2Et] V2Et]-2U4Et] V4Et])~4 (-mU2Et] 

V2pE0]+m U4Et] V4pE0])~2 (((-U2pE0] V2Et]+U4pE0] V4Et]) (-m U2Et] V2pE0]+m 
U4Et] V4pE0]))/(2 U2Et] V2Et]-2 U4Et] V4Et])-2-(m (U2p EO] U4Et]-U2Et] U4pE0]) 
(V2pE0] V4Et]-V2Et] V4pE0]))/(2 U2 Et] V2Et]-2 U4Et] V4Et])''2)) +(4 m vep2i 
veqli (-U2pEt] U4Et]+U2Et] U4pEt])~2 (U2 Et] (V2p EO] V3Et]-V2Et] V3p EO])+U1 Et] 
(-V2pE0] V4Et]+V2Et] V4pE0])+U4Et] (V3pE0] V4Et]-V3Et] V4pE0])))/((2 U2 Et] 

V2Et]-2 U4Et] V4Et]) (-m U2Et] V2p EO]+m U4Et] V4p EO] ) "2) - (8 m sqlqli (-U2pEt] 
U4Et]+U2Et] U4pEt]) (UlpEt] U2Et] V2Et]-UlEt] U2pEt] V2Et]+U2pEt] U4Et] V3Et]-U2Et] 
U4pEt] V3Et]-UlpEt] U4Et] V4Et]+uiEt] U4pEt] V4Et]) (U2Et] (V2pE0] V3Et]-V2Et] 
V3pE0])+UlEt] (-V2pE0] V4Et]+V2Et] V4pE0])+U4Et] (V3pE0] V4Et]-V3Et] V4pE0])))/((2 
U2Et] V2Et]-2 U4Et] V4Et])''2 (-m U2 Et] V2p EO]+m U4Et] V4pE0]))-(8 m veqli~2 
(-U2pEt] U4Et]+U2Et] U4pEt]) (UlpEt] U2Et] V2Et]-uiEt] U2pEt] V2Et]+U2pEt] 

U4Et] V3Et]-U2Et] U4pEt] V3Et]-UlpEt] U4Et] V4Et]+uiEt] U4pEt] V4Et]) (U2Et] 

(V2pE0] V3Et]-V2Et] V3pE0])+UlEt] (-V2pE0] V4Et]+V2Et] V4pE0])+U4Et] (V3pE0] 
V4Et]-V3Et] V4pE0])))/((2 U2Et] V2 Et]-2 U4Et] V4Et])''2 (-m U2Et] V2p EO]+ni 
U4Et] V4pE0]))-(8 m veqli veq2i (-U2pEt] U4Et]+U2Et] U4pEt]) (U2pEt] U4Et] 
VlEt]-U2Et] U4pEt] VlEt]-U2pEt] U3Et] V2Et]+U2Et] U3pEt] V2Et]-U3pEt] U4Et] 
V4Et]+U3Et] U4pEt] V4Et]) (U2Et] (V2pE0] V3Et]-V2Et] V3pE0])+UlEt] (-V2pE0] 
V4Et]+V2Et] V4pE0])+U4Et] (V3pE0] V4Et]-V3Et] V4pE0])))/((2 U2Et] V2 Et]-2 
U4Et] V4Et])''2 (-m U2Et] V2p EO]+m U4Et] V4pE0])) + (4 m sqlpli (-U2pEt] U4Et]+U2Et] 
U4pEt]) (U2pEt] V2Et]-U4pEt] V4Et]) (U2Et] (V2pE0] V3Et]-V2Et] V3pE0])+UlEt] 
(-V2pE0] V4Et]+V2Et] V4pE0])+U4Et] (V3p EO] V4Et]-V3Et] V4pE0])))/((2 U2 Et] 
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V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(4 m vepli 
veqli (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2 [t] (V2p[0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~3 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) 
+ (8 m sqlqli (U2p [0] U4[t]-U2[t] U4p [0]) (U2p [t] V2[t]-U4p[t] V4[t]) (Ulp[t] 

U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] 

U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (U2[t] (V2p [0] V3[t]-V2[t] V3p [0])+U1 [t] 

(-V2p[0] V4[t]+V2[t] V4p[0])+U4[t] (V3p [0] V4[t]-V3[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m 
U4[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p[0]) 
(V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2)) +(8 m veqli~2 
(U2p[0] U4[t]-U2[t] U4p[0]) (U2p[t] V2[t]-U4p[t] V4[t]) (Ulp[t] U2[t] V2[t]-Ul[t] 
U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] 

U4p[t] V4[t]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0] )+Ul [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2[t]-2 U4[t] V4[t])~4 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(8 m veqli veq2i (U2p [0] U4[t]-U2[t] 
U4p[0]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] 
U3[t] V2[t]+U2[t] U3p[t] V2[t]-U3p[t] U4[t] V4[t]+U3[t] U4p[t] V4[t]) (U2[t] 

(V2p[0] V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p[0] 
V4[t]-V3[t] V4p[0])))/((2 U2[t] V2 [t]-2 U4[t] V4[t])-4 (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+ni U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4 [t] ) "2-(m 
(U2p[0] U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])~2)) +(8 m sqlqli (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p [t] V2[t]-U4p[t] 
V4[t]) (U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] 
(U2p[0] V2[t]-U4p[0] V4[t])) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0] )+Ul [t] (-V2p[0] 
V4[t]+V2[t] V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2 [t]-2 
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U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(8 m veqli~2 (-U2p[t] U4[t]+U2[t] 
U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] 
(-U2p[0] V3[t]+Ulp[0] V4[t])+Ul[t] (U2p[0] V2[t]-U4p[0] V4[t])) (U2 [t] (V2p [0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) 

+ (3 m veqli veq2i (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) 

(U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] V4[t])+U3[t] 

(U2p[0] V2[t]-U4p[0] V4[t])) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0] )+Ul [t] (-V2p[0] 
V4[t]+V2[t] V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2 [t]-2 
U4[t] V4[t])~4 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(3 m vep2i veqli (U2p [0] U4[t]-U2[t] 
U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) (U2[t] (V2p [0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~3 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(3 m''2 vep2i veqli (U2p [0] U4[t]-U2[t] 
U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t])-2 (V2p [0] V4[t]-V2[t] V4p [0]) (U2[t] 

(V2p[0] V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] 
V4[t]-V3[t] V4p[0])))/((2 U2[t] V2 [t]-2 U4[t] V4[t])~3 (-m U2[t] V2p [0]+m 
U4[t] V4p[0])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(4 m''2 sqlpli (-U2p[t] U4[t]+U2[t] 
U4p[t])~2 (V2p[0] V4[t]-V2[t] V4p [0] ) (U2 [t] (V2p [0] V3[t]-V2[t] V3p [0])+U1 [t] 
(-V2p[0] V4[t]+V2[t] V4p[0])+U4[t] (V3p [0] V4[t]-V3[t] V4p [0]) ) ) / ( (2 U2 [t] 

V2[t]-2 U4[t] V4[t])~3 (-m U2 [t] V2p[0]+m U4[t] V4p [0]) (((-U2p[0] V2[t]+U4p[0] 
V4[t]) (-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])''2-(m 
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(U2p[0] U4[t]-U2[t] U4p[0]) (V2p [0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] V2[t]-2 
U4[t] V4[t])''2)) -(4 m~2 vepli veqli (-U2p[t] U4[t]+U2[t] U4p[t])''2 (V2p[0] 
V4[t]-V2[t] V4p[0]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2[t]-2 U4[t] V4[t])~3 
(-m U2[t] V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) -(8 
m~2 sqlqli (U2p [0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p [t]) (Ulp[t] 

U2[t] V2[t]-Ul[t] U2p[t] V2[t]+U2p[t] U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] 

U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (V2p[0] V4[t]-V2[t] V4p [0]) (U2 [t] (V2p [0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(8 m''2 veqli''2 (U2p [0] U4[t]-U2[t] 
U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (Ulp[t] U2 [t] V2[t]-Ul[t] U2p [t] V2[t]+U2p[t] 
U4[t] V3[t]-U2[t] U4p[t] V3[t]-Ulp[t] U4[t] V4[t]+Ul[t] U4p[t] V4[t]) (V2p [0] 
V4[t]-V2[t] V4p[0]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2[t]-2 U4[t] V4[t])~4 
(-m U2[t] V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2 [t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) -(8 
m~2 veqli veq2i (U2p [0] U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) 

(U2p[t] U4[t] Vl[t]-U2[t] U4p[t] Vl[t]-U2p[t] U3[t] V2[t]+U2[t] U3p [t] V2[t]-U3p[t] 
U4[t] V4[t]+U3[t] U4p[t] V4[t]) (V2p[0] V4[t]-V2[t] V4p [0]) (U2 [t] (V2p [0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0])))/((2 U2[t] V2[t]-2 U4[t] V4[t])~4 (-m U2 [t] V2p [0]+m U4[t] V4p [0]) 
(((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2 [t] 

V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])-2)) -(8 ni~2 sqlqli (-U2p[t] U4[t]+U2[t] 
U4p[t])~2 (U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] 
V4[t])+Ul[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p [0] V4[t]-V2[t] V4p [0]) (U2[t] 
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(V2p[0] V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] 
V4[t]-V3[t] V4p[0])))/((2 U2[t] V2 [t]-2 U4[t] V4[t])''4 (-m U2[t] V2p [0]+ni 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(8 m~2 veqli''2 (-U2p[t] U4[t]+U2[t] 
U4p[t])~2 (U2[t] (-Ulp[0] V2[t]+U4p[0] V3[t])+U4[t] (-U2p[0] V3[t]+Ulp[0] 

V4[t])+Ul[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p [0] V4[t]-V2[t] V4p [0]) (U2[t] 

(V2p[0] V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] 
V4[t]-V3[t] V4p[0])))/((2 U2[t] V2 [t]-2 U4[t] V4[t])''4 (-m U2[t] V2p [0]+ni 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) -(8 ni''2 veqli veq2i (-U2p[t] U4[t]+U2[t] 
U4p[t])~2 (U2[t] (U4p[0] Vl[t]-U3p[0] V2[t])+U4[t] (-U2p[0] Vl[t]+U3p[0] 

V4[t])+U3[t] (U2p[0] V2[t]-U4p[0] V4[t])) (V2p[0] V4[t]-V2[t] V4p [0]) (U2[t] 

(V2p[0] V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] 
V4[t]-V3[t] V4p[0])))/((2 U2[t] V2 [t]-2 U4[t] V4[t])''4 (-m U2[t] V2p [0]+ni 
U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] V2p [0]+m U4 [t] V4p[0]))/(2 
U2[t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] 
V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])''2)) +(8 ni''2 veqli veq2i (-U2p[t] U4[t]+U2[t] 
U4p[t])~2 (U2[t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] 
(-V2p[0] V4[t]+V2[t] V4p[0])) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0] )+Ul [t] (-V2p[0] 
V4[t]+V2[t] V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2 [t]-2 
U4[t] V4[t])''2 (-m U2[t] V2p[0]+m U4[t] V4p [0])''2)-(16 m~2 veqli veq2i (U2p [0] 
U4[t]-U2[t] U4p[0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] V2[t]-U4p[t] V4[t]) 

(U2[t] (-Vlp[0] V2[t]+Vl[t] V2p[0])+U4[t] (Vlp[0] V4[t]-Vl[t] V4p [0] )+U3 [t] 

(-V2p[0] V4[t]+V2[t] V4p[0])) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0] )+Ul [t] (-V2p[0] 
V4[t]+V2[t] V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 U2[t] V2 [t]-2 
U4[t] V4[t])~4 (-m U2[t] V2p [0]+m U4 [t] V4p [0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) 

(-m U2[t] V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] 
U4[t]-U2[t] U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) 

+ (16 m~3 veqli veq2i (U2p [0] U4[t]-U2[t] U4p [0]) (-U2p[t] U4[t]+U2[t] U4p[t])''2 
(V2p[0] V4[t]-V2[t] V4p[0]) (U2 [t] (-Vlp[0] V2[t]+Vl[t] V2p [0] )+U4 [t] (Vlp[0] 
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V4[t]-Vl[t] V4p[0])+U3[t] (-V2p[0] V4[t]+V2[t] V4p [0])) (U2[t] (V2p [0] V3[t]-V2[t] 
V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0] )+U4 [t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ) / ( (2 
U2[t] V2[t]-2 U4[t] V4[t])''4 (-mU2[t] V2p [0]+m U4 [t] V4p[0])''2 (((-U2p[0] 
V2[t]+U4p[0] V4[t]) (-m U2 [t] V2p [0]+m U4[t] V4p[0]))/(2 U2[t] V2 [t]-2 U4[t] 
V4[t])"2-(m (U2p[0] U4[t]-U2[t] U4p [0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2 [t] 
V2[t]-2 U4[t] V4[t])''2)) +(4 m~2 sqlqli (-U2p[t] U4[t]+U2[t] U4p[t])~2 (U2 [t] 

(V2p[0] V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] 
V4[t]-V3[t] V4p[0]))~2)/((2 U2[t] V2[t]-2 U4[t] V4[t])~2 (-m U2 [t] V2p [0]+m 
U4[t] V4p[0] )''2) + (4 m~2 veqli~2 (-U2p[t] U4[t]+U2[t] U4p[t])''2 (U2 [t] (V2p [0] 
V3[t]-V2[t] V3p[0])+Ul[t] (-V2p[0] V4[t]+V2[t] V4p [0])+U4 [t] (V3p [0] V4[t]-V3[t] 
V4p[0]))~2)/((2 U2[t] V2[t]-2 U4[t] V4[t])~2 (-m U2 [t] V2p[0]+m U4[t] V4p [0]) ~2) - (8 
m~2 sqlqli (U2p [0] U4[t]-U2[t] U4p [0]) (-U2p[t] U4[t]+U2[t] U4p [t]) (U2p[t] 
V2[t]-U4p[t] V4[t]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ~2) / ( (2 U2 [t] V2[t]-2 U4[t] V4[t])~4 
(-m U2[t] V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) -(8 
m~2 veqli~2 (U2p [0] U4[t]-U2[t] U4p [0]) (-U2p[t] U4[t]+U2[t] U4p[t]) (U2p[t] 
V2[t]-U4p[t] V4[t]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ~2) / ( (2 U2 [t] V2[t]-2 U4[t] V4[t])~4 
(-m U2[t] V2p[0]+m U4[t] V4p[0]) (((-U2p[0] V2[t]+U4p[0] V4[t]) (-m U2 [t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2[t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) +(8 
ni~3 sqlqli (U2p [0] U4[t]-U2[t] U4p [0]) (-U2p[t] U4[t]+U2[t] U4p[t])~2 (V2p [0] 
V4[t]-V2[t] V4p[0]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ~2) / ( (2 U2 [t] V2[t]-2 U4[t] V4[t])~4 
(-m U2[t] V2p[0]+m U4[t] V4p[0])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] 

V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2 [t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)) +(8 
ni~3 veqli~2 (U2p[0] U4[t]-U2[t] U4p [0]) (-U2p[t] U4[t]+U2[t] U4p[t])~2 (V2p[0] 
V4[t]-V2[t] V4p[0]) (U2[t] (V2p[0] V3[t]-V2[t] V3p [0])+U1 [t] (-V2p[0] V4[t]+V2[t] 
V4p[0])+U4[t] (V3p[0] V4[t]-V3[t] V4p [0] ) ) ~2) / ( (2 U2 [t] V2[t]-2 U4[t] V4[t])~4 
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(-m U2[t] V2p[0]+m U4[t] V4p[0])~2 (((-U2p[0] V2[t]+U4p[0] V4[t]) (-mU2[t] 
V2p[0]+m U4[t] V4p[0]))/(2 U2 [t] V2 [t]-2 U4[t] V4[t])~2-(m (U2p [0] U4[t]-U2[t] 
U4p[0]) (V2p[0] V4[t]-V2[t] V4p[0]))/(2 U2[t] V2[t]-2 U4[t] V4[t])~2)); 
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